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ABSTRACT 

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of 
their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the 
radio. X-ray, and infrared (IR) wavebands. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts 
from the AGN and Galaxy Evolution Survey (AGES). We select AGNs with observations in the radio at 1 .4 GHz 
from the Westerbork Synthesis Radio Telescope, X-rays from the Chandra XBootes Survey, and mid-IR from 
the Spitzer IRAC Shallow Survey. The radio. X-ray, and IR AGN samples show modest overlap, indicating 
that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy 
colors and luminosities, as well as Eddington ratios, for obscured or optically faint AGNs. We also measure the 
two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 Mpc, and derive typical dark 
matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red sequence galaxies, are 
strongly clustered (with Mhaio ^ 3 x lO'-' Mq), and have very low Eddington ratios (A < lO"-'); (2) X-ray- 
selected AGNs are preferentially found in galaxies that lie in the "green valley" of color-magnitude space and 
are clustered similar to typical AGES galaxies (Mhdo ^ lO'-' Mq), with 10"^ < A < 1; (3) IR AGNs reside 
in slightly bluer, slightly less luminous galaxies than X-ray AGNs, are weakly clustered (Mhaio 10'^ /;"' 
Mq), and have A > 10"^. We interpret these results in terms of a simple model of AGN and galaxy evolution, 
whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches 
a critical mass between ~10'^ and lO'"* Mq. After this event, star formation ceases and AGN accretion shifts 
from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically faint, radio-bright) modes. 

Subject headings: galaxies: active — quasars: general — large-scale structure of universe — radio continuum: 
galaxies — surveys — X-rays: galaxies 



1. INTRODUCTION 

There is increasing evidence that the evolution of galax- 
ies is related to the evolution of their supermassive black 
holes (SMBHs), which primarily grow through accretion of 
material as active galactic nuclei (AGNs). This connec- 
tion between galaxies and SMBHs is suggested by the ob- 
served tight c orrelation betwee n SMBH and galaxy bulge 
masses (e.g., [M agon'ia n et al.l Il998t iFerrarese & Merritj 
l2000HGebhardt et al. 200df, and the similar redshift evolution 
of star form ation and AGN activity (e.g., 'Madau et al. "1996|; 
lUeda et an[2003 ). Galaxies show a well-established bimodal- 
ity in color that separates blue, generally disk-dominated, 
star-forming galaxies (the "blue cloud") from red, gener- 
ally bulge-dominated galaxies (the "red sequence") (see e.g., 
IStrateva eta l."2001: Blanton 2006; Faber et al. 2007, and ref- 
erences therein). Recent work has suggested that AGNs may 
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play a crucial role in the origin of this bimodality, and par- 
ticularly the quenching of star formatio n in blue galaxies and 
their transition to the red sequence (e.g..|Hopkins et al.l2006al : 
ICroton et al.l[2006l: iKhalatyan et allllOOSlP 

Valuable clues to the associated evolution of galaxies 
and AGNs come from the properties of AGN hosts. It is 
well-established that radio AGNs are m ost commonly found 
in massive earl y-type galaxies (e.g., lYee & GreenI 119871 : 



iBest et al] l2005h . and that a majority of bright local elhp- 
ticals show low-power nuclear radio sources (Sadler et ajj 
Il989l). Among optically selected AGNs, Kauffmann et al] 
(12003a) found that low-redshift, low-luminosity AGNs 
reside in massive galaxies with generally old stellar pop- 
ulations, and that the stellar population becomes younger 
for galaxies with increasing [O JII] luminosity. At higher 
redshifts (0.5 < z < 1.5), X-ray-selected AGNs are also 
found in massive galaxie s with generally red colors 
(lAlonso-Herrero et al.l l2008l). although some studies have 
also found an excess of AGNs in the "green valley" between 
the blue cloud and red sequence, compared to the distribution 
JSanchez et al. 



of qu i escent galaxies 



2004; NandraetalJ 
2007; Georgakakis et al. 2008; Silverman et al. 200a 
Kocevski et al. 2009b; Schawinski et al. 2009b; Treister et al] 
2009). The differences in host galaxies between AGNs of 
different types suggest that the process of nuclear accretion is 
related to the mass and age of the host galaxy. 

Further constraints on the links between AGNs and galax- 
ies come from the clustering of AGNs and their local galaxy 
environments. Among quiescent galaxies, red galaxies are 
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more strongly clustered, and therefore are found in denser 
environm ents, than blue ga l axies (for some recent results 
see e.g., IZehavi et al.l 120051; ICoil et al.1 l2008h . Clustering 
measurements on scales > 1 Mpc allow us to estimate the 
masses of the dark ma tter halos in which galaxies reside (e.g., 
ISheth & TormenI [l999 ). There is evidence that the evolu- 
tion of galaxies is strongly linked to host halo mass (e.g., 
iFaber et al.l2007h . and similar relationships may also exist for 
AGN activity. Previous works have found somewhat differing 
results for AGN clustering and env ironments, la rgely owing 
to different samples of sources (see lBrown et aLl uOOl, for a 
summary of previous results). However, a general picture has 
emerged in which radio-loud AGNs are found in dense envi- 
ronments, and massive halos, similar to galaxy clusters, while 
radio-quiet AGNs are found in poorer environments similar to 
field galaxies. 

Recently, wit h large surveys such as the Sloan Digital Sky 
Surv ev (SPSS. lYork et al.l l200Q). the 2dF QSO R edshift Sur- 
vey ( Croom et al .1 12004 . and the DEEP2 survey dPavis et alJ 
l2003h . AGN-galaxy clustering results have been placed on 
firmer statistical footing. Studies of the autocorrelation of 
optically selected quasars from 2dF and SDSS showed that 
their clustering amplitude increases with redshift out to z '--^ 3, 
indicating that at all redshifts they inhabit dar k matter halos 
of similar characteristic mass ~3 x 10'^ (Porciani et al." 
2004; Croom et al. 2005; Coil et al. 200 7 ; Myers et al. 2007,; 
Shen et allllOOTt Ida Angela etll] 120081; iPadmanabhan et all 
2008h . This suggests that luminous quasar activity is related 
to the mass of the surrounding dark matter halo. Among 
lower luminosity optically selected AGNs from SDSS, the 
Seyfert galaxi es are generally less c lustered than the full 
galaxy sample (iKauffmann et al.lf2004l) . but show comparable 
large-scale clustering to control samples of q uiescent galax- 
ies with similar propertie s to the AGN hosts jLi et al.l 120061; 
iMandelbaum et al.' 2009). Very low luminosity AGNs (low 
ionization nuclear emission regions (LIN ERs)) show no sig- 
nificant bias relative to no rmal galaxies^ jMiller et al.ll2003|; 
IConstantin & Vogelevll2006i) . At higher redshifts. ICoil et al I 
(I2007h found that quasars in the DEEP2 fields at z ~ 1 show 
a slight antibias (although at the l-2cr level) relative to all 
galaxies, similar to the antibias observed for blue galaxies. 

In contrast, using X-ray da ta in the Extended Grot h Strip 
region and DEEP2 redshifts, Georgakakis et al.l (l2007h found 
that X-ray selected AGNs generally lie in overdense regions 
at z ^ 1, w hile in the Extended Chandra Deep Field South, 
ISilvermari et al. (2008) found that AGNs in "green" hosts at 
0.4 < z < 1.1 preferentially reside in large-scale structures 
such as walls or filaments. 

X-ray AGNs are generally strongly clustered and reside 
in relatively dense environments. Studies of the spatial 
auto-correlation of X-ray AGNs detected with Chandra and 
XMM-Newton (e.g., iB asilakos et al. 2004; Gilli et al. 20051 
lYang et al.l2006HPuccetti et al.,.2006:Miyaji et al...200 7) gen- 
erally obtain relatively large clustering amplitudes (although 
with differing results that may reflect a variation in cluster- 
ing with X-ray flux; Plionis et al. 2008). Recently, Coil et al. 
(|2009) found that X-ray AGNs in the Extended Groth Strip 
are significantly biased relative to galaxies in the DEEP2 sur- 
vey, i n contrast to the results fo r optical quasars. Correspond- 
ingly, iGeorgakakis et al.l (12007 ) showed that X-ray AGNs in 
DEEP2 at z 1 generally lie in overdens e regions, while 
in the Extended Chandra Deep Field South. [Silverman et al.l 
(l2008h found that AGNs in "green" hosts at 0.4 < z < 



1.1 preferentially reside in large-scale structures such as 
walls or filaments. O t her recent works (iLehmeret al] 120091; 
lKocevskietaLll2009al; iGalametz etakl l2009l) have found an 
enhancement of the density of X-ray selected AGNs in large- 
scale structures at relatively high redshifts (z > 0.9). How- 
ever, the significance of the results for X-ray AGNs has been 
limited by relatively small AGN samples. 

Compared to X-ray selected sources, radio-selected AGN 
are found in even dens er environments i ncluding X-ray 
groups and clusters (e.g., Croston et al. 2008). Radio AGNs 
show strong spat ial clustering similar to local ellipti cal galax- 
ies (e.g. Mande lbaum et al.ll2009l: IWake et alJl2008l) . with ev- 
idence for higher am plitude for more powerful sources (e.g., 
lOverzier et al1l2003h . 

These results indicate that AGNs selected using different 
techniques represent separate populations, with differences in 
host galaxies, environments, and accretion modes. Large red- 
shift surveys with extensive multiwavelength coverage now 
make it possible to explore, within the same uniform data set, 
(1) host galaxies, (2) clustering, and thus dark matter halo 
masses, and (3) Eddington ratios of AGNs selected by these 
different techniques. 

In this paper, we study a sample of AGNs in the red- 
shift interval 0.25 < z < 0.8, using spectroscopic redshifts 
from the AGN and Galaxy Evolution Survey (AGES). We fo- 
cus on AGNs selected in the radio. X-ray, and infrared (IR) 
bands, which provide relatively efficient and unbiased ways 
of iden tifying certain populations of AGNs (for a review see 
iMusho tzky 2004). We do not, in general, use optical spectral 
diagnostics for selecting AGNs, because of selection effects 
that can arise from the low-resolution fiber spectroscopy of 
some of the sources. We place particular emphasis on ob- 
scured or optically faint AGNs, for which we can study the 
optical properties of their host galaxies. For these objects we 
can also measure host galaxy bulge luminosities and thus esti- 
mate black hole masses (assuming an Lbui-A^BH relation), in 
order to derive Eddington ratios. 

This paper is organized as follows. In Section |2] we de- 
scribe the Bootes multiwavelength data set. In Section [3] we 
describe the main AGES galaxy sample, and in Section |4] we 
discuss the samples of radio. X-ray, and IR selected AGNs. 
In Section|5]we present results on the colors and luminosities 
of AGN host galaxies. We describe our correlation analysis 
in Section |6] and present the results, along with estimates of 
dark matter halo mass, in Section [T] In Section |8] we calcu- 
late average X-ray spectra and Eddington ratios for the dif- 
ferent samples of AGNs. In Section |9] we discuss our results 
in terms of a simple picture of AGN and galaxy evolution, 
in which luminous quasar activity and the termination of star 
formation occur at a characteristic dark matter halo mass. We 
summarize our results in Section[TOl We also include two ap- 
pendices that describe (Appendix lAb our correction for AGN 
contamination of host galaxy colors, and (AppendixlBl) limits 
on the contamination of the IR AGN sample by star-forming 
galaxies. 

Throughout this paper we assume a lambda cold dark mat- 
ter (ACDM) cosmology with fin, = 0.3 and ^lf^ = 0.7. For 
direct comparison with other works, we assume Hq = 70 km 
s"' Mpc"', except for comoving distances and absolute mag- 
nitudes, which are explicitly given in terms of /i = //o/( 1 00 km 
s-i Mpc"'). The NOAO Deep Wide-Field Survey (NDWFS) 

' Throughout this paper, the terms "quiescent" or "normal" refer to galax- 
ies with no detected AGN activity, but do include star-forming galaxies. 
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photometry in the Bw, R, and / bands is presented in Vega 
magnitudes, while magnitudes in the SDSS u and r bands are 
AB normalized. All quoted uncertainties are la (68% confi- 
dence). 

2. OBSERVATIONS 

The 9 deg^ survey region in Bootes covered by the ND- 
WFS (iJannuzi & Devi figg g*) is unique among extragalactic 
multiwavelength surveys in its wide field and uniform cov- 
erage using space- and ground-based observatories, including 
the Chandra X-Ray Observatory and the Spitzer Space Tele- 
scope. Extensive optical spectroscopy makes this field espe- 
cially well suited for studying the statistical properties of a 
large number of AGNs (C. Kochanek et al. 2009, in prepara- 
tion). The areas covered by the X-ray, optical, infrared, and 
radio surveys are shown in Figure 1 . 

Redshifts for this study come from AGES, which 
used the Hectospec multifiber spectrograph on the MMT 
dFabricant et al.ll2005l) . We use AGES Data Release 2 (DR 2), 
which consists of all the AGES spectra taken in 2004-2006. 
The 2004 observations covered 15 pointings for a total area 
of 7.9 deg^, which we define to be the main AGES survey 
region. In the 2004 observations, targets include (1) all ex- 
tended sources (i.e. galaxies, see Section |4~4| | with R < 19.2 
(2) a randomly selected sample of 20% of all extended sources 
with 19 < < 20, (3) all extended sources with R <20 and 
Infrared Arary Camera (IRAC) 3.6, 4.5, 5.8, and 8.0 ^im 
magnitudes <15.2, 15.2, 14.7, and 13.2, respectively. In 
addition, (4) fainter sources were observed, selected mainly 
from objects that are counterparts of Chandra X-ray sources 
dMurray et al. 2005; Brand et al. 2006a; Kenter et al. 2005), 
radi o sources from the Very Large Array (VLA) FIRST sur- 
vey ( iBecker et al.ll 19951) . and objects selected from 24 fim ob- 
servations with the Multiband Imaging Photometer for Spitzer 
(MIPS) (E. Le Floc'h et al. 2009, in prepai-ation). AGES DR 
2 contains /-selected targets with / < 21.5 for point sources 
and / < 20.5 for extended sources. 

AGES redshifts and spectral classifications follow the pro- 
cedures used by SDSS. Redshifts are measured by comparing 
the extracted spectroscopy to the grid of galaxy templates, and 
choosing the combination of template and redshift that mini- 
mizes the between the data and the model. Good redshifts 
are defined to be those that do not have a second minimum 
in that is within 0.01 in reduced x^ to the best-fit value. 
All redshifts were examined by eye and bad fits were flagged 
manually. 

Optical photometry from NDWFS was used for the selec- 
tion of AGES targets and to derive optical colors and fluxes 
for AGES sources. NDWFS images were obtained with the 
Mosaic- 1 camera on the 4-m Mayall Telescope at Kitt Peak 
National Observatory, with 50% completeness limits of 26.7, 
25.0, and 24.9 mag, in the Bw, R, and / ba nds, respectively. 
Photo metry is derived using SExtractor dBertin & ArnoutsI 
fT996h . 

Radio data are taken from the Westerbork Synthesis Radio 
Telescope (WSRT) 1.4 GHz radio survey that covers ^ 7 deg^ 
in the center of the NDWFS field ( de Vries et all 2002). 3035 
radio sources are detected in the AGES survey region, to a 
limiting flux of «0.1 mJy and beam size 13" x 27". Centroid 
positions of the radio sources are measured to 0"4. 

X-ray data are taken from the XBootes survey, which cov- 
ers the full AGES spectroscopic region. The XBootes mo- 
saic consists of 126 5 ks Chandra ACIS-I exposures and is 
the largest contiguous field observed to date with Chandra 
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Fig . 1 . — Map of the Bootes survey region, showing the approximate areas 
covered by the AGES, WSRT (radio), XBootes (X-ray), and IRAC Shallow 
Survey (mid-IR) observations. 

(iMurray et al.ll2005h . Owing to the shallow exposures and 
low background in the ACIS detector. X-ray sources can be 
detected to high significance with as few as four counts. The 
XBootes catalo g contains 3 293 X-ray point sources with four 
or more counts CKenter et al . 2005), and of these, 2724 lie in 
the AGES survey region. 

Mid-IR ob servations are taken from the Spitzer IRAC Shal- 
low Survey (lEisenhardt et al.l l2004l) . which covers the full 
AGES field in all four IRAC bands (3.6, 4.5, 5.8, and 8 /im), 
with 5ct flux limits of 6.4, 8.8, 51 and 50 /iJ y respectively. 
The IR AC photometry is described in detail in lBrodwin et al.l 
( |2006|) . The catalog contains 15,488 sources detected in all 
four IRAC bands, and of these, 14,069 lie within the AGES 
survey region. 

3. GALAXY SAMPLE 

The first step in our analysis is to define a statistical sample 
of galaxies from the AGES catalog. These galaxies will be 
used to provide a comparison sample against which we can 
compare the properties of AGN host galaxies. These galax- 
ies will also be used to determine the spatial cross-correlation 
between AGNs and normal galaxies, and so to derive AGN 
clustering. 

The statistical galaxy sample consists of AGES galaxy tar- 
gets with the flux limit / < 20. We outline the key aspects 
of the statistical galaxy sample here, while full details will 
be presented by D. Eisenstein et al. (2009; in preparation). 
AGES galaxy targets must satisfy quality cuts in the / band 
and either the Bw and R bands in the NDWFS imaging. In 
addition, all objects in the sample must have good detections 
(SExti-actor FLAGS < 8) in all three bands, to allow for K 
corrections (this eliminates very few sources). In addition, ob- 
jects must be detected as an extended source (see Section l4~4] i 
in at least one of the three bands. 

In this analysis we use the auto (Kron) /-band magni- 
tude; however in some regions the /-band photometry suffers 
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Fig. 2. — Positions of 6262 AGES galaxies in the main statistical sample 
(Section |3) at redshifts 0.25 < z < 0.8. Points are color-coded by redshift. 
The black line shows the boundaiy of the 15 main AGES pointings. 

from low-surface-brightness halos around bright stars, which 
causes the Kron galaxy magnitudes to be overly bright. To 
correct for this, we define the Ir magnitude constructed from 
the sum of the R band Kron magnitude and the R-I color de- 
rived from 6" aperture magnitudes. We then compare Ir and 
the / (Kron) magnitude and compute /tot, which consists of 
the fainter magnitude (if the two differ significantly) and oth- 
erwise is the average of the two magnitudes, /tot is typically 
0.01-0.09 mag fainter than the / (Kron) magnitude. We cut 
the sample at /tot < 19.95, which excludes 2% of the galaxies 
with Kron magnitude / < 20. The positions on the sky of the 
AGES galaxies with 0.25 < z < 0.8 are shown in Figure |2] 
and their redshift distribution is shown in Figure [3] 

3.1. Sampling weights 

For the targets that pass the above criteria, a complicated set 
of random sparse sampling criteria were used to select objects 
for observation. In general, 100% of bright (/ < 18.5) galax- 
ies were observed, while the sampling rate for sources with 
18.5 </< 20 was 20%-30% (depending on whether they 
pass flux cuts in interesting UV, optical, or IR bands). In ad- 
dition, many galaxies were assigned fibers (and obtained red- 
shifts) as part of the AGES observations, but were not system- 
atically targeted. Some of these redshifts correspond to X-ray 
or radio-selected AGNs and so are used in this analysis. How- 
ever, these objects are not included in the "main" statistical 
galaxy sample that we use for AGN-galaxy cross-correlation. 
The main AGES sample includes only those galaxies that 
were systematically targeted in the random sampling. 

For each target in the main AGES sample that has a well- 
determined redshift, we assign a series of sampling weights. 
First we assign a sparse sampling weight, equal to the inverse 
of the random sampling rate that is known exactly from the 
target selection (that is, if 20% of objects in a class are tar- 
geted, the sampling weight is 5). Second, we determine a tar- 
get assignment weight, given by the inverse of the probability 
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Fig. 3. — Redshift distribution for galaxi es in the main AGES sample and 
for the AGES random catalog (see Section IT4l . The solid and dotted black 
histograms show the redshift distribution for AGES main sample galaxies (in 
our analysis we only include objects at 0.25 < z < 0.8, shown by the solid 
black histogram; for completeness the dotted black histogram shows AGES 
main sample galaxies at z < 0.25). The red dashed histogram shows the 
redshift distribution for the random galaxy catalog. For clarity, the histogram 
values for the random galaxies are divided by 200, since the random sample 
contains 200 times as many objects as the (statistically weighted) observed 
galaxy sample. Note that the observed galaxy histogram is not normalized by 
statistical weights. 



that a target was assigned a fiber in the 2004-2006 observing 
runs. Owing to the high completeness of the AGES observing 
strategy for the main galaxy sample, this fiber weight aver- 
ages only 1.05. Finally, we assign a redshift weight, given 
by the inverse probability that an object that is assigned to a 
fiber succeeds in yielding a well-determined redshift. We de- 
termine this weight for various subsamples of the targets in 
bins of magnitude, surface brightness, color, and sparse sam- 
pling class, by computing the ratio of the number of attempts 
to the number of well-determined redshifts. The total statis- 
tical weight W associated with each galaxy is the product of 
these three weights. 

3.2. Bright star exclusion 

Owing to difficulties in obtaining redshifts for objects in 
the vicinity of bright stars, objects around bright sources were 
not targeted for AGES spectroscopy. The exclusion was per- 
formed using a bright source list taken from objects in the 
USNO-B catalog (Monet et al. 2003) with R<n. A number 
of these bright "stars" correspond to extended sources, some 
of which are AGES target galaxies. To minimize the exclu- 
sion of targets around AGES galaxies, the 475 objects that 
match extended sources in the Two Micron All Sky Survey 
(2MASS) catalog (Skrutskie et al. 2006) were removed from 
the bright star list before the AGES targets were selected. This 
leaves 6690 bright sources around which objects were not tar- 
geted for AGES spectroscopy. Possible target objects were 
excluded around these bright sources, with an exclusion ra- 
dius that varied with the flux of the bright source and the flux 
of the target object. Around a bright source of magnitude R, 
faint target objects were excluded out to a radius 

ninft = 20"+5"(i5-/;), (1) 

while brighter target objects were excluded inside smaller 
radii. 
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To minimize biases in our spatial clustering analysis, we 
avoid the regions on the sky from which AGES targets were 
excluded. For simplicity, we remove from our sample all 
AGES sources (regardless of flux) within ru^ii of a bright 
source, including the bright source itself if it is an AGES 
galaxy. We exclude the same regions in our random cata- 
log (Section [3.4l i. This exclusion removes 776 AGES target 
galaxies that correspond to objects in the bright source cata- 
log, but were not matched to 2MASS extended sources. How- 
ever, all but 19 of these bright AGES targets lie at z < 0.25. To 
minimize the effects of removing the bright AGES galaxies, 
we restrict our analysis to objects at z > 0.25, which leaves 
us with 6262 galaxies, with total statistical weight equal to 
15,653. We refer to these galaxies as the "main" AGES sta- 
tistical sample and focus on this sample in our analysis. 

3.3. K-corrections and colors 

One powerful diagnostic for understanding the properties 
of galaxies is their distribution in color and absolute mag- 
nitude. Many previous studies have found that the color- 
magnitude distribution of galaxies is bimodal, separated into 
a red sequence of luminous, bulge-dominated, passively 
evolving galaxies, and a blue cloud of less luminous, star- 
forming, disk-dominated galaxies (e.g.. IStrateva et al.ll2001l: 
iBlanton et a l. 2003b': iKauffmann et al.ll2003d) . 

For the AGES galaxies, we use the 4" aperture magnitudes 
and /tot to perform /T-corrections and thus derive the rest- 
frame colors and absolute magnitudes. As discussed above, 
only galaxies with good Bw, R, and / photometry are included 
in the galaxy sample. We obtain Bw and R magnitudes by 
adding the Bw-I and /?-/ 4" aperture colors, respectively, to 
/tot. These magnitudes are in general close to the Kron mag- 
nitudes; the mean and standard deviation in Am = mtot-»!Kron 
are (0.01, 0.10) and (-0.003, 0.038) in the Bw and R bands, 
respectively. We convert Vega magnitudes to AB using cor- 
rections of H-0.02 mag, H-0.19 mag, and H-0.44 mag for the Bw, 
R, and / bands, respectiv ely, and correct absolute magnitudes 
for interstellar extinction (iSchlegel et al.ll99&) . which for this 
field is only A/ = 0.02 mag. 

We define rest-frame colors and absolute magnitudes in 
terms of photometric bands that are close to the wavelengths 
probed by the NDWFS photometry at typical AGES redshifts. 
We convert the / and Bw magnitudes to ° 'r and ° 'm respec- 
tively, which are defined as the AB magnitudes in the SDSS 
r and u bands, shifted blueward by z = 0.1. These conver- 
sions are insensitive to the galaxy spectral energy distribution 
(SED) at z = 0.42 and 0.27, respectively. Calculating rest- 
frame colors in terms of these bands allows f or straightfor- 
ward c om parison to results from S PSS (e.g.. iBlanton et al.l 
I2003allbll5 IKauffmann et aT]|2003bl) . The jjT-corr ections are 
performed using kcorrect version 3.2 (IBlanton et al. 2003a); 
using an updated version of kcorrect (version 4. 1 .4) produces 
no significant change in the /^-corrections. 

In addition to calculating rest-frame colors for each galaxy, 
we derive the absolute magnitude (Moi,.) in the ° V band at 
z = 0.1. We also calculate the absolute magnitude corrected 
for passive evolution, assuming 1.6 mag of evolution per 
unit r edshift relative to the observed redshift (lEisenstein et al.l 
l2005l) . This correction varies from at z = 0.1 (by definition) 
to H- 1 . 1 mag at z = 0.8, in the sense that the evolution-corrected 
absolute magnitude at z = 0. 1 is fainter than the observed Moa^. 
The evolution-corrected absolute magnitudes are used only in 
the estimates of galaxy bulge and black hole masses (Sec- 
tion [82|i; elsewhere we use absolute magnitudes that are not 



corrected for evolution. 

Figure |4] shows rest-frame " '(w-r) color versus Moi^ for 
galaxies in the AGES main sample in the redshift interval 
0.25 < z < 0.8. The bimodal distribution in galaxy color (the 
red sequence and blue cloud) is clearly evident. We divide 
the galaxies into red and blue subsamples; to account for the 
slope of the red sequence, we define the quantity 

A = °\u-r)+Om(MoA, + 20), (2) 

(where Mo i,. is calculated for h= I) and divide red and blue 
galaxies by a cut in A. The red sequence evolves slightly in 
color with redshift, so we empirically determine the median 
value in A for the red sequence as a function of z; this varies 
from Ared = 2.64 at z = 0.25 to Ared = 2.52 at z = 0.8. At each 
redshift, we set the boundary between red and blue samples 
to be 0.3 mag blueward of A,ed. The selection boundaries for 
z = 0.25 and z = 0.8 are shown in Figured This criterion se- 
lects 3119 red galaxies and 3143 blue galaxies in the redshift 
interval 0.25 < z < 0.8. After a small correction for nuclear 
contamination from AGNs (AppendixlAli. we select 3146 red 
galaxies and 3119 blue galaxies. 

3.4. Random galaxy catalog 

To derive the spatial correlation of galaxies and AGNs (Sec- 
tion|6]l requires an estimate of the number of galaxy-galaxy or 
AGN-galaxy pairs expected if the galaxies were randomly dis- 
tributed in comoving space. We therefore produce a catalog 
of galaxies that are spatially uncorrected, but reflect the se- 
lection function of the AGES sample. To ensure that the frac- 
tional uncertainties of the AGN-random catalog pairs are neg- 
ligible, the density of random galaxies must be significantly 
larger than in the data catalog. This is especially true on small 
scales (less than 1 /i"' Mpc), where real galaxies are most 
strongly clustered. We design the random catalog to contain 
200 times the (statistically weighted) number of galaxies in 
the data catalog. 

For each galaxy in the AGES statistical sample with sta- 
tistical weight W, we generate A^,- duplicate galaxies, where 
A^, = 200W, rounded to the nearest integer. We then assign 
each of these random galaxies a statistical weight equal to 
Wjf = W /Ni, so that the total statistical weight in the random 
sample is equal to the total weight in the data. 

We place each random galaxy in a random sky position in- 
side the main AGES field (and outside the bright star masks) 
and offset its redshift by a random -0.05 < Az < +0.05. This 
effectively applies a boxcar smoothing to the line of sight dis- 
tances on scales of 200-300 Mpc, larger than the largest 
observed structures, and so eliminates any significant corre- 
lation along the line of sight. The redshift distribution of the 
random catalog is shown in Figure [3] along with the distribu- 
tion for the observed galaxies. 

4. AGN SAMPLE 

We next define the samples of AGNs for which we will de- 
termine host galaxy, clustering, and accretion properties. For 
complete as possible an AGN census, we select AGNs inde- 
pendently in the radio. X-ray, and IR wavebands. As dis- 
cussed above, we do not use optical spectroscopic selection 
in order to avoid systematic effects in the analysis of the fiber 
spectroscopy. We restrict ourselves to sources that have red- 
shifts from AGES in the interval 0.25 < z < 0.8, are located 
in the main AGES field, and are away from bright stars, as de- 
scribed in Section [3] In this section we describe the selection 
criteria and AGN samples. 
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Fig. 4. — (a) Rest-frame optical luminosities and colors for AGES galaxies. The contours and points show the distribution of rest-frame (u — r) color versus 
Mo.ii- absolute magnitude (not corrected for passive evolution), for AGES galaxies in the statistical sample in the redshift interval 0.25 < z < 0.8 (points are 
left out in the region covered by contours). The diagonal dashed lines show the empirical division between red and blue galaxy samples for z = 0.25 (top) and 
z = 0.8 (bottom), (b) Distiibution in the color parameter A (see Eqn.|2) for AGES galaxies at 0.25 < z < 0.8, showing the "red sequence" and "blue cloud". The 
boundaries between red and blue galaxies at z = 0.25 and z = 0.8 are shown as dashed lines, as in (a). The dotted line shows the distribution including galaxies at 
z < 0.25, which include many lower luminosity blue galaxies, and which more clearly shows the "green valley". 
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Fig. 5.— Positions of 585 AGNs in the main AGES field. We show AGNs 
selected in the radio (orange circles) X-rays (green stars), and IR (red squares) 
that lie in the main AGES field and away from bright stars, and in the redshift 
interval 0.25 < z < 0.8. The solid line shows the boundary of the main AGES 
field. 



For each radio source we calculate the radio power f i 4 ghz 
(in W Hz"'). To obtain the rest-frame power we include a 
small /T-correction of 0. 1 dex, appropriate for the average red- 
shift of the radio-detected objects between 0.25 < z < 0.8, 
and the typical spec trum of faint 1 .4 G Hz sources (a w 0.5, 
where Si, oc Pra ndoni et al.ll206 6h. The distributions in 
A .4 GHz and redshift are shown in Figure |6}j. For objects 
with low radio power, the observed emission can be pow- 
ered by either AGN activity or star formation. It is possi- 
ble to separate AGNs from star-forming galaxies based on a 
combination of Pi.4 ghz and optical spectral features such as 
Ha luminosity (e.g. jKauffmann et al.ll2008l) : however the re- 
quired fits to the optical spectra for the AGES galaxies are 
not currently available. To minimize the number of star- 
forming galaxies, we restrict our sample to the 122 sources 
with log Pi .4 GHz > 23.8. In low -redshift samples, almost all 
galaxies with radio power above this cut are AGNs (see Fig- 
ure 1 of iKauffmann et alj |2008|) '°. All the radio AGNs in 
our sample have Pi .4ghz < 2 x 10^^ W Hz~', and so a re in 
the luminosity range typical of Fan aroff & RilevI d 1 9741) type 
I sources (Ledlow & Owen 1996). The positions on the sky 
of the radio AGNs are shown in Figure |5] To test for spuri- 
ous matches of the radio sources to AGES galaxies, we offset 
the radio source positions by 1' and reperformed the source 
matching. We find only four matches, for which we calcu- 
lated Pi .4 GHz corresponding to the redshift of the matched 
AGES galaxy. None of the matches have logPi 4ghz > 23.8, 
indicating that there is minimal contamination from spurious 
matches in our radio AGN sample. 



4.1. Radio AGNs 

We first select AGNs in the radio. The AGES catalog con- 
tains 525 matches (within 2") to WSRT radio sources, of 
which 245 are in the AGES field and in the interval 0.25 < 
z < 0.8 (note that the WSRT field is «7 deg^ and does not 
cover the full AGES spectroscopic region, as shown in Fig- 
ure [T]). None of the radio sources match more than one AGES 
galaxy. 



4.2. X-ray AGNs 

We next define the X-ray AGN sample. Among the 3293 X- 
ray sources with >4 X-ray counts in the iKenter et alJ (l2005h 
catalog, 362 of these lie within the AGES field, are not close 
to bright stars (see Section [372] |. and are matched within 3. "5 
to objects with good AGES redshifts at 0.25 < z < 0.8. We 



Note that IKauffmann et alj 120081) assume Hq = 
cut is scaled to our assumed cosmology. 



50 km s Mpc ; our 
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have excluded four objects that show evidence for extended 
X-ray emission that could arise from galaxy groups or clusters 
(we verify explicitly that this does not significantly affect the 
clustering results). Nine of the 362 sources in our X-ray AGN 
sample have two AGES galaxies within 3 "5; for these we se- 
lected the closer galaxy as the optical counterpart. We test 
for spurious matches by offsetting the IR-selected AGN po- 
sitions by 1' and reperforming the source matching. We find 
five matches, indicating contamination from spurious matches 
of only wl%. 

For each X-ray source we calculate an 0.5-7 keV luminos- 
ity from the number of detected source counts in the 0.5-2 
keV and 2-7 keV bands separately, and assuming a power- 
law X- ray spectrum with photon index F = 1.8, typical of 
AGNs (To zzi et akl 12006). The corresponding conversions 
from count rates (in counts s"') to flux (in erg cm"^ s~') are 
5.9 X 10"'^ erg cm"^ count"' in the 0.5-2 keV ba nd and 1 .9 x 
10~" erg cm"^ count"' in the 2-7 keV band ( Renter et alj 
l2005h . For simplicity we do not include /^-corrections, which 
are generally small (<5%); for most sources uncertainty in 
the flux and F are significantly larger. 

The X-ray AGN positions are shown in Figure |5] and their 
distribution in luminosity and redshift is shown in Figure |6l7. 
The luminosities of these sources are in the range 10^^ <Lx < 
10^^ erg s"', characteristic of moderate- to high luminosity 
AGNs and significantl y larger than is typic ally found for star- 
forming galaxies (e.g- lRanalli et al.ll2003h . We are therefore 
confident that the X-ray emission from these objects is pow- 
ered by nuclear accretion rather than star formation. 

We note that since XBootes sources have as few as 4 counts 
in 5 ks, the fluxes of the faintest sources are affected by 
the Eddington bias. For a differential flux distribution with 
dN/dS oc S~^, the true mean flux for a source obse rved with 
n counts is approximately n — f3 (lEddingtonl 1 1 9 1 3l) . At the 
XBootes flux limits (5o.5 -2keV ^ 5 x 10" '^ erg s"'), /3 ~ 2 
dVikhlinin et al. 1995; Ke nteret al"1l2005h . so a so urce with 
4 cou nts has a true mean flux of ~2 counts (iKenter et al.l 
l2005h . This correction is relatively small for most sources, 
and does not affect our interpretation of the X-ray sources as 
accretion-dominated. Correcting for the Eddington bias, all 
of the sources have Lx > 10'*^ erg s~', so we conclude that 
contamination by starburst-dominated galaxies will be small. 
In what follows we use Lx values that are corrected for Ed- 
dington bias. 

4.3. IRAGNs 

We finally select AGNs using IRAC observations. With 
the launch of Spitzer, the mid-IR has provided a sen- 
sitive new window for identifying AGNs. The mid-IR 
SEDs of broad-line quasars and Seyfert galaxies typically 
show a featureless, roughly power-law continuum that rises 
to longer wavelengths (e.g., iGlikman et al.l 120061) . This 
SED can be used to separate AGNs from normal and star- 
burst galaxies, whose SEDs exhibit a stellar bump at 1.6 
fim and then fall at longer wavelengths. Several crite- 
ria have been dev eloped to select AGNs based on IRAC 
color s (iLacv et al. 2004; Stern et al. 2005) or SED fitting 
(e.g., lAlonso-Herrero et al. 2006' Donlev et al.. .2007). I n 
this paper we use the color selection of IStern et al.l (l2005h . 
which was derived using the Bootes multiwavelength data 
and so is well matched to the depth of the survey. Be- 
cause mid-IR wavelengths are not as strongly affected by 
dust extinction as the optical or UV, IRAC selection is par- 
ticularly useful for identifying dust-obscured AGNs (e.g.. 
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Fig. 6. — Distribution in redsliift and luminosity for (a) radio, (h) X-ray, 
and (c) IR AGNs witli AGES redsliifts, in tlie main AGES field and away from 
bright stars. Luminosities are the observed 1.4 GHz radio power (Pi.4gHz): 
the 0.5-7 keV luminosity (Lx), and uL,y at 4.5 fj,m n the rest frame, (L4.5^m)- 
We focus on the redshift interval 0.25 < z < 0.8, shown by the vertical dot- 
ted lines, and in the radio, sources with logPi 4ghz > 23.8 (see Section|4), 
shown by a horizontal line. Objects included in the analyses are shown by 
filled symbols. The dashed lines show approximate luminosity Hmits for the 
survey flux limits (note that the X-ray limit is for an X-ray spectrum with 
photon index F = 1.8). The histograms show the distribution in z and lumi- 
nosity for sources in this redshift interval. For the radio sources, the dashed 
histograms show the distributions in z and Pi 4 qhz f™' the complete sample 
of radio sources in the interval 0.25 < z < 0.8, while the solid histograms 
show those sources with logPi 4ghz > 23.8. 
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Stern et al.ll2005l: iRowan-Robinson et ani2005l: [Poiletta et alj 
2006HHickox et alJ[ 2007 [) . 

The IStern et alJ ( 2005 ) AGN color-color criteria requires 
that a source be detected in all four IRAC bands. Of the IRAC 
sources detected in all four bands, 2952 are matched (within 
3"5) to AGES sources at 0.25 < z < 0.8. The distribution 
in IRAC colors for these objects are shown in Figure |7] The 
objects in the bottom left of this diagram have IRAC colors 
similar to stars ([5.8] -[8.0] < 0.8 and [3.6]- [4.5] < 0.3); 
these are generally passively evolving galaxies that are dom- 
inated by starlight in the mid-IR. The diagram also shows a 
continuous sequence of sources with redder [5.8] - [8.0] and 
[3.6] -[4.5] colors, which are generally star-forming galax- 
ies with dust features (mainly due to polycyclic aromatic hy- 
drocarbons (PAHs)) that move through the IRAC bands with 
redshift (see Figure 1 of Stern et al. 2005, for evolutionary 
tracks of star-forming galaxies on this diagram). A separate 
sequence of sources with very red [3.6]- [4.5] colors is se- 
lected by the IStern et al.l (l2005h AGN criteria (shown by the 
dashed line). These SEDs are inconsistent with typical emis- 
sion from dust heated by star formation, and are characteristic 
of broad-line quasars and Seyfert galaxies. 

Of the AGES sources at 0.25 < z < 0.8 matched to the 
IRAC photometry, these criteria select 238 objects. Six of 
these sources have two AGES galaxies within 3. "5; we se- 
lected the closer galaxy as the optical counterpart. We test 
for spurious matches by offsetting the IR-selected AGN po- 
sitions by 1' and find six matches, indicating contamination 
from spurious matches of only w2%. The positions of the IR- 
selected AGNs are shown in Figure |5] For each IR-selected 
AGN, we calculate the rest-frame luminosity (i^Lij) at 4.5 //m, 
L4.5fj.m- The SED at rest-frame 4.5 fim is probed by the IRAC 
photometry at all redshifts in our sample, and since this wave- 
length is significantly redward of the 1 .6 /im peak in the stel- 
lar emission from galaxies, the observed flux is likely domi- 
nated by the AGN. We estimate -L4,5pm by logarithmically in- 
terpolating between the observed lyL^, in the 4.5 /im, 5.8 /im, 
and 8 /im IRAC bands. The distribution of the IRAC-selected 
AGNs in ^4.5^™ and redshift is shown in Figure|6]:. 

One possible limitation of IRAC color selection is in reli- 
ability. For the radio and X-ray samples, the high radio and 
X-ray luminosities are robust indicators of the presence of nu- 
clear accretion. The rehability of IRAC selecti on is less well 
studie d [alt hough there ar e disc ussions of this in lBarmby et alJ 
(|2006l) and iGoriian et all (120081) ]. The most likely contami- 
nants are star-forming galaxies, which in some extreme cases, 
can have IR luminosities and colors similar to our IR-selected 
AGN sample. 

Our sample of AGNs with AGES spectroscopy appears to 
suffer little contamination from star-forming galaxies, owing 
prim arily to the r elatively bright optical and IRAC flux lim- 
its. iDonlev et al.l (12008) examined Spitzer 24 /im sources 
in the GOODS-S field, and found that IRAC color-selected 
samples in deep surveys have significant contamination from 
star-forming galaxies, but th is contaminat i on de creases for 
higher IR and optical fluxes. iDonlev et alj (l2008h argue that 
the selection of AGNs by fitting a power-law SED to the 
observed IRAC emission reduces the contamination in their 
sample. To check whether power-law selection would im- 
prove the rehabihty of our sample, we apply the power-law 
criteria of iDonley et al.l (|2008) to objects with AGES spectra 
at 0.25 < z < 0.8 and detections in all four IRAC bands. We 
identify 54 power -law AGNs , of w hich all but one are also 
selected using the lStern et alJ (l2005h color selection. Among 



1.0 



be 
> 



0.5 



0.0 



-0.5 



o Radio 
♦ X-ray 




6 X 



-0.5 0.0 0.5 1.0 1.5 2.0 
[5.8]-[8.0] (Vega) 



2.5 



3.0 



Fig . 7 . — IRAC color-color diagram for all AGES sources at 0.25 < z < 0.8 
that lie inside the main AGES field and away from bright stars, and have 5cr 
detections in all four IRAC bands. The dashed line shows the selection region 
defined by Stern et al. (2005). Symbols for radio and X-ray AGNs are as in 
Figure [3] Most radio AGNs are found to the lower left, in the region occu- 
pied by stars and quiescent galaxies I Stern et al. 2005), while X-ray AGNs 
are found thi'oughout the distribution, with a significant number in the IRAC 
AGN selection region. 

the power-law AGNs, 53 lie in the X-ray field, and 21 (40%) 
are detected in X-rays, indicating that they are unambiguously 
AGNs. In comparison, 184 color-selected AGNs are not se- 
lected by the power-law criteria and have X-ray coverage, and 
of these 97 (53%) are detected in X-rays. Therefore in our 
sample, the IRAC color criteria are equally or more likely to 
identify X-ray-detected AGNs as the power-law criteria. We 
conclude that for our relatively shallow IRAC o bservations 
and br ight optical spectroscopic flux limit, the IStern et al.l 
(120051) color criteria are as reliable as power-law criteria in 
selecting AGNs. 

As another check, we carefully analyze the average X-ray 
emission from the IRAC color-selected sources that are not 
detected in X-rays. The average fluxes, stacked spectrum, 
and distribution of observed X-ray counts indicate that star- 
forming contamination in our IR-selected AGN sample is at 
most 20%, and most likely <10% (see Appendix IB1 for a de- 
tailed discussion). We conclude that contamination of the 
IRAC AGN sample from star-forming galaxies is small. 

4.4. Optical counterparts 

Useful observational clues as to the nature of AGNs come 
from the morphologies of their optical counterparts. Objects 
that are unresolved in the NDWFS images are generally dom- 
inated in the optical by the active nucleus, while objects that 
are optically extended are dominated by their host galaxies 
and have optically faint or obscured nuclei. For optically ex- 
tended sources, we can determine the optical properties of the 
host galaxies, and indirectly, estimate their central black hole 
masses. We therefore divide the AGN sample into two cat- 
egories based on whether their optical counterparts are ex- 
tended or unresolved. 

The sample of AGNs with extended counterparts is defined 
as sources having the SExtractor parameter CLASS_STAR < 
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0.8 in all of the Bw, R, and / bands. 92%, 97%, and 99% of 
these objects have CLASS_STAR < 0.3 in the Bw. R, and / 
bands, respectively. These will be referred to as AGNs with 
"galaxy" (or "extended") counterparts. Almost all (113) of 
the radio AGNs, roughly two-thirds (238) of the X-ray AGNs, 
and half (133) of the IR AGNs have extended optical counter- 
parts. For simplicity, and to obtain the best possible statistics, 
we include all the optically extended sources in our analy- 
sis. We note that, 74, 157, and 111, of the radio. X-ray, and 
IR AGNs, respectively were systematically targeted by AGES 
and are included in the main galaxy sample. We verify explic- 
itly that if we limit the sample only to these sources (which 
may suffer less complicated selection effects), there is no sig- 
nificant change to our results. 

Although the optical counterparts of these sources are ex- 
tended, the derivation of host galaxy properties is complicated 
by the fact that an optically weak AGN can still affect the in- 
tegrated colors from the galaxy. By comparing the aperture 
photometry for these AGNs to that of normal galaxies, we 
carefully estimate the AGN color contamination and derive 
corrections to the host galaxy colors. Details of these cor- 
rections are given in Appendix lAl The color corrections are 
small (typically ~0.1 mag), and do not significantly affect our 
conclusions. 

The second category of AGN optical counterparts are those 
with pointlike optical morphologies. We define pointlike 
sources as those with CLASS_STAR > 0.8 in any of the Bw, 
R, and / bands. The flux from most of these sources is dom- 
inated by nuclear emission in the Bw', for simplicity, we will 
refer to them as "unresolved", but we note that many of these 
sources are extended in the R and / bands in which there is 
a larger contribution from the host galaxy. Although these 
objects are dominated by the nucleus at blue optical wave- 
lengths, they would not historically be defined as "quasars". 
Their absolute magnitudes in the B band (estimated roughly 
from their observed 7? -band flux) are - 1 8 < Mb < -24. Above 
the usual luminosity threshold defined for quasars (Mb < 
-23), the unresolved AGN sample only includes 1 1 sources 
(or 23 sources if we include unresolved objects not selected 
in the radio. X-ray, or IR). Therefore, our sample of optically 
unresolved AGNs would primarily be classified as luminous 
Seyfert galaxies. We note that it is not surprising that there are 
few luminous optical quasars in our sample at 0.25 < z < 0.8. 
From t he optical quasar luminosity function of iRichards et al.l 
(l2005h . we would expect only 18 objects with Mb < -23 in 
this redshift interval over the area of 7.3 deg- (after excluding 
bright stars), which probes a volume of 5 x 10^ Mpc^ at 
0.25 <z< 0.8. 

4.5. Overlap between samples 

One complication of selecting AGNs in different wave- 
bands is that the samples overlap. Some AGNs are selected 
in more than one waveband, so the different AGN samples 
are not entirely independent. The various samples and their 
overlaps are given in Table [1] and a Venn diagram depicting 
the relative sizes and overlaps of the samples is shown in Fig- 
ure |8] This figure shows the sample of sources that lie in the 
area covered by the radio. X-ray, and IRAC observations, to 
allow a fair comparison between the samples. 

It is immediately apparent that there is very little overlap 
between the radio AGNs and those selected in the X-ray and 
IR. This may provide an important clue as to the nature of the 
accretion in these systems (Section|9]l. However, there is sig- 
nificant overlap between the X-ray and IR samples. Roughly 



one-third of X-ray AGNs are in the IR sample, and half of the 
IR AGNs are detected in X-rays. These fractions are smaller 
(roughly 20% and 40%, respectively) when we consider only 
optically extended AGNs. The sources that are selected as 
both X-ray and IRAC AGNs have relatively high luminosi- 
ties; their median Lx and ^4.5^™ are ~50% higher than for 
AGNs selected in only the X-ray and IR bands, respectively. 
At some stages in the following analysis we study the proper- 
ties of this joint X-ray and IR sample, but in most of the paper 
we consider the X-ray and IR AGN samples separately. This 
allows our results to be compared directly to other works that 
study X-ray or IR AGNs only, but it is important to keep in 
mind that the samples are not entirely independent. 

It is also essential to stress that AGN selection is strongly 
affected by the flux limits of the observations in each wave- 
band. In all three bands, our AGN selection criteria identify 
objects that are dominatedhy the AGN in that band. The lumi- 
nosity cut in the radio, the X-ray flux limit, and the color cri- 
teria in the IR all select objects where the nucleus dominates 
over stellar processes at the observed wavelengths. However, 
as discussed in detail by Ho (2008), if observations go deep 
enough, then the majority of AGNs show signatures of nu- 
clear activity in multiple bands, particularly in the X-ray and 
radio. Indeed, we show in AppendixiBjthat many, if not most, 
of IR AGNs that are undetected in the shallow X-ray survey 
still have low-level X-ray emission. 

We also note that most X-ray AGNs are not selected by 
the IRAC color-color criteria, but in these objects the IRAC 
emission from the host galaxy may mask a relatively low- 
luminosity, red power-law-type SED from the AGN. Figure|7] 
shows that in the IRAC color-color diagram, many X-ray 
AGNs lie just below the IStern et aP (i2b05b selection region, 
which could be evid ence for a low-level AGN contribution 
jGoriian et al.l l2008h . It is possible that these objects con- 
tain relatively IR-faint AGNs that redden the colors of a qui- 
escent galaxy and move it toward (but not all the way into) 
the AGN color-color s election reg i on. Ho wever, IRAC color- 
color space below the IStern et all ( l2005b region is also pop- 
ulated by quiescent galaxies with "green" colors similar to 
many X-ray AGN hosts (Section |5]l, so the colors alone do 
not provide evidence for a weak AGN. In the future, a detailed 
study of the multiwavelength SEDs of these sources could put 
limits on the nuclear contribution to the IRAC flux. 

5. AGN HOST GALAXIES 

Links between AGN activity and the evolution of galaxies 
are reflected in the characteristics of the galaxies that host dif- 
ferent classes of AGNs. As discussed in Section |373l a power- 
ful diagnostic for understanding galaxy properties is their po- 
sitions in optical color-magnitude space. In this section, we 
examine the optical colors and luminosities for AGNs with 
extended optical counterparts (corrected for nuclear contami- 
nation as described in Appendix IaIi. and compare to the dis- 
tribution for the full galaxy sample. 

The color-magnitude diagrams for AGN hosts are shown 
in Figure |9}j. AGNs of all types have optical colors and lu- 
minosities that lie within, or close to, the parameter space 
occupied by normal galaxies, however their distributions in 
color-magnitude space are markedly different from the total 
galaxy population, and vary significantly for AGNs selected 
in different wavebands. 

Radio AGNs are concentrated on the luminous end of the 
red sequence, indicating that they reside in massive early- 
type galaxies, with only a few sources detected in the blue 
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TABLE 1 
AGN SAMPLES 









Also selected in 






Sample 


Total" 


Only'' Radio= 


X-ray"^ 


IRAC^^ 


AH 






All AGNs 










Radio 


122 (122) 


103 (103) 


6 (6) 


6 


(6) 


7 


X-ray 


362 (296) 


238 (192) 6 (6) 




111 


(91) 


7 


IRAC 


238 (199) 


114 (95) 6 (6) 


111 fPiJ 






7 






Optically extended AGNs 








Radio 


113 (113) 


98 (98) 


4 r^j 


5 


(5) 


6 


X-ray 


241 (200) 


181(153) 4 (4) 




44 


(37) 


6 


IRAC 


133 (114) 


78 (66) 5 f5j 


44 (37) 






6 






Optically unresolved AGNs 








Radio 
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Note. — Values in parentheses are for sources inside the region covered by radio. X-ray, 
and IRAC observations (Fig.lTi. 

^ Total number of AGNs with AGES redshifts 0.25 <z< 0.8, selected in waveband for that 
row, inside the main AGES field, away from bright stars.'' The number of AGNs selected in the 
waveband for that row only.^ The number of AGNs selected in the waveband for that row, and 
also in the band for that column only.'^ The number of AGNs selected in all three wavebands. 




Fig. 8. — Venn diagram showing the relative number of AGNs with AGES 
redshifts at 0.25 < z < 0.8, selected in the radio, X-ray, and IR (see also 
Table [T). The figure shows only those sources that lie within the regions 
covered by radio. X-ray, and IRAC observations (Figure[T). The overlapping 
areas between the samples correspond to the relative numbers selected in 
multiple wavebands, except for the region where all three overlap, which is 
slightly too small in the diagram. There is some overlap between the X-ray 
and IRAC samples, but little overlap between these and the radio AGNs. 



cloud (Figure |9] top panels). The mean and dispersion in 
° '(M-r) and Mo.i, are ((C), tic) = (2.6,0.2) and ((M),(Tm) = 
(-22.0,0.6), respectively. In contrast. X-ray AGNs are found 
throughout the galaxy color-magnitude space, with a dispro- 
portionate number found in the "green valley" (Figure|9] cen- 
ter panels). X-ray AGN hosts have ((C), tic) = (2.2,0.4) and 
((M),crM) = (-21.3,0.7). The fraction of AGES main sample 
galaxies containing X-ray AGNs is 4.6% ±0.7% in the "green 
valley" (having A within ±0.2 mag of the division between 
red and blue galaxies). In contrast, the X-ray AGN fraction 
is 2.1% ± 0.2% for galaxies with redder or bluer colors. As 
we discuss below, this enhanced X-ray AGN activity in green 
galaxies may reflect the connection between these objects and 
SMBH accretion. 

IR-selected AGNs have a color distribution similar to that 
of X-ray AGNs, although they have a less pronounced excess 
in the "green valley", and have very few hosts on the red se- 
quence (Figure|9] bottom panels), with {{C),ac) = (2.1,0.4). 
The hosts of IR AGNs have smaller average optical luminosi- 



ties than those X-ray AGNs, with {{M),(7m) = (-21.1,0.6). 
The sources that are selected as both X-ray and IRAC AGNs 
are found in a similar region of color-magnitude space to the 
full X-ray and IR selected samples, with ((C), tic) = (2.2,0.4) 
and ((M),crM) = (-21.2,0.6). 

We conclude that IR and radio AGNs occupy almost en- 
tirely separate regions of the color-magnitude space, with ra- 
dio AGNs in luminous red sequence galaxies, and IR AGNs 
in the less luminous blue and "green" hosts. X-ray AGNs are 
found in galaxies that are slightly redder and more luminous 
than the hosts of IR-selected AGNs. 

6. CORRELATION ANALYSIS 

Another important clue to the evolution of AGNs and their 
host galaxies comes from the large-scale environments in 
which AGNs are found. We derive the environments of AGNs 
and galaxies through a spatial two-point correlation analysis, 
which probes the clustering of sources on different comoving 
scales. We first calculate the autocorrelation of AGES main 
sample galaxies (and the red and blue samples separately), 
to determine the absolute linear bias of AGES galaxies rela- 
tive to dark matter. Using the absolute bias, we then derive 
the characteristic masses of the dark matter halos which host 
galaxies of different types. 

We next derive the cross-correlation of AGNs with AGES 
galaxies, to determine the relative bias between AGNs and 
galaxies, from which we derive the absolute bias and char- 
acteristic dark matter halo masses for the AGNs themselves. 
AGN-galaxy cross-correlation has been used since the first 
studies of quasar clustering (e.g., Bahcall et al. 1969) and has 
two main advantages, as discussed in lCoil et al.l (I2007D . First, 
it does not require a full understanding of the AGN selection 
function, which may be complicated, particularly for AGNs 
selected in multiple wavebands. Instead, we need only know 
the selection function for the galaxies, which is well deter- 
mined (Section[3]). The second advantage is that AGN-galaxy 
cross-correlation has much greater statistical power than an 
AGN autocorrelation, owing to the large number of galax- 
ies in the redshift survey (greater than 6000) compared to the 
number of AGNs (a few hundred). 

The two-point correlation function ^(r) is defined as the 
probability above Poisson of finding a galaxy in a volume el- 
ement dV at a physical separation r from another randomly 
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Fig. 9. — (a) Optical colors and absolute magnitudes of AGNs with extended optical counterparts. Contours, black points, and dashed Hnes are as in Figure|4] 
Orange circles, green stars, and red squares show radio. X-ray, and IR AGNs, respectively. The arrows show the typical correction for nuclear contamination 
for different host galaxy colors, which range from to 0.3 mag (Appendix|A). (h) Distribution in A for AGNs selected in the three wavebands, compared to all 
AGES main sample galaxies at 0.25 < z < 0.8 (scaled by 1/25, thick gray line). The radio AGN color distribution peaks along the red sequence, while X-ray 
AGNs are found preferentially in the "green valley" between the red sequence and the blue cloud. The distribution of IR AGNs is similar to that of X-ray AGNs, 
although they are typically found in somewhat less luminous galaxies with very few on the red sequence, and show a less pronounced peak in the "green valley". 



chosen galaxy, such that 

dP = n[l + ^(r)]dV, (3) 

where n is the mean space density of the galaxies in the sam- 
ple. To calculate th e autocorrelation of AG ES galaxies, we 
derive ^(r) using the lLandv & Szalavl(ll993h estimator: 



where DD, DR, and RR are the number of data-data, data- 
random, and random-random galaxy pairs, respectively, at a 
separation r. We account for the sparse sampling of fainter 
galaxies by multiplying each pair by the product of the sta- 
tistical weights of the observed (or random) galaxies. That 
is, when correlating galaxy sample 1 against galaxy sample 2, 



ar)=-^(DD-2DR+RR), 



(4) 
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the weighted number of data-data pairs is 



DD: 



iGDiD: 



wiwi, 



(5) 



where Wl and W2 for each pair are the statistical weights of 
the galaxies from samples 1 and 2, respectively. By including 
these weights, we ensure that brightest galaxies do not overly 
dominate the correlation function. We c orre spondingly in- 
clude the random galaxy weights (Section |331l in calculating 
DR and RR. The random galaxy weights are assigned so that 
their total weight equals that of the AGES galaxies, so that 
DD, DR, and RR may be comp ared directly. 

Following C oil et all (l2007b . for the cross-correlation be- 
tween AGNs and galaxies we use the simple estimator 

D1D2 



D1R2 



■1 



(6) 



where D1D2 is the number of AGN-galaxy pairs, and D1R2 
is the number of AGN-random pairs. Because the selection 
function of the AGN samples is not as well defined as that 
for the AGES main sample galaxies, we do not assign statis- 
tical weights to the AGNs. Therefore we use only the galaxy 
weights in calculating the weighted number of pairs, such that 



D1D2 



!GDiD2 



(7) 



and correspondingly for D1R2. 

In the range of separations I ^ r < IQ /i"' Mpc, ^(r) 
for galaxies is roughly observed to be a power law, ^(r) = 
{r/r())~^ , although recent work has shown evidence for sep- 
arate terms in the correlation owing to dark matter halos 
that host single galaxies and those that host pairs of galaxies 
(e.g., Zehavi et al."2004; 'Zheng et al."2007|; ICoil et al.ll2008t 
|Brown et al. 2008; Zheng et al. 2009). For simplicity, and in 
light of the uncertainties in the AGN-galaxy correlations, we 
will restrict our fits to power-law models. Throughout the pa- 
per ro is given in comoving coordinates. 

6.1. The projected correlation function 

With a redshift survey, we cannot directly measure ^(r) in 
physical space, because peculiar motions of galaxies distort 
the line-of-sight distances measured from redshifts. To ac- 
count for these redshift-space distortions, we instead measure 
the correlation function ^(rp,7r), where r,, and it are the pro- 
jected comoving separations between galaxies in the direc- 
tions perpendicular and parallel, respectively, to the mean line 
of sight between the two galaxies. 

Redshift space distortions only affect the correlation func- 
tion along the line of sight, so we integrate £,{rp,Tr) along the 
77 direction to obtain the statistic. 



Wpirp) = 2 / (,(rp,n)dn, 
Jo 



(8) 



whi ch is independent of re dshift space distortions (follow- 
ing |D^ir&Peeblei[l983|). This estimator has been used 
to measure co rrelations in a number of survey s, for ex- 
ample SDSS (iZehavi et all llOOl iLi et all l2006l). 2 SLAO 
jWake et al.' '2008*). DEEP2 jCoil et al.1 l2007l |2008"1). and 
GOODS (Gilli et al. 2007). If 

'''max = 00, then we average over 
all line-of-sight peculiar velocities, and Wp{rp) can be directly 
related to ^(r) (for a power-law parametrization) by 



Wp(rp) = rp 



r(l/2)r[(7-l)/2] 

r(7/2) 



(9) 



In practice, we truncate the integral at a finite TTmax, to re- 
duce the noise from integrating over weakly correlated objects 
at large tt. Here we use Tr^ax = 25 Mpc. Because of the 
finite TTmax, in order to recover ro and 7 exactly we must ac- 
count for redshift-space distortions due to peculiar motions 
along the line of sight. This requires an accounting for coher- 
ent infall of galaxies on large scales, as well as the random 
velocities of galaxies on small scales ("fingers of God"). 

We acco unt for th ese red shift-s pace distortions using the 
method of Coil et all ('2007, 2008). Coherent infall depends 
on the parameter /? « /b (Ka^iser 1987), where ilm evalu- 
ated at the mean redshift of the sample, and b is the absolute 
bias relative to dark matter. For z = 0.5 (the mean redshift for 
the AGNs), = 0.59. We take b = 1.23, the absolute bias 
of AGES galaxies at 0.25 < z < 0.8 (see Section|63]and Sec- 
tionlTJ, so that (3 = 0.59. For the pairwise random velocities of 
galaxies, we take ai2 = 500 km s"', which is characteristic of 
galaxies with luminosities typical of the AGES sample (e.g., 
iHawkins et alll2003t ICoil et alJl2008h . For reasonable varia- 
tions in /3 (±0.1) or an (±400 km s"'), the resulting Wp(rp) 
changes by at most a few percent. 

We then derive ^(rp,7r) in redshift space as a function of 
ro and 7, including coher ent infall and random velocities as 
described in Section 4. 1 of IHawkins et al.l (l2003h . We numer- 
ically integrate £,{rp,TT) to TTmax = 25 Mpc to determine a 
model Wp(rp) which we fit to the data. For a given ro and 7, 
this Wpifp) differs from that given in Eqn.|9]by less than 4% 
for rp < 10 Mpc. Accordingly, the best-fit ro and 7 differ 
from those using Eqn.|9]by less than a few percent. 

6.2. Calculating relative bias 

For each subset of galaxies and AGNs, we also calculate the 
relative bias brei compared to the autocorrelation of all AGES 
main sample galaxies. The relative bias for the red galaxy 
sample is given by 

, w,,(rr,) [red -red] 

wp(r^) [gal-gal]' ^'^^ 

and accordingly for blue galaxies. For AGNs, since we calcu- 
late the AGN-galaxy cross-correlation, the relative bias is 

Wp(r^)[AGN-gal] 

Orel — ■ ^) 

Wp(rp) [gal -gal] 
In each case we derive b^ei as a function of rp, and calculate 
the mean biases averaged over scales of 1-10 Mpc and 
0.3-1 /!-' Mpc. 

The galaxy autocorrelation varies with redshift, owing to 
the evolution of large-scale structure, and because the use of a 
flux-limited sample means we selected more luminous galax- 
ies at higher z- This can affect bias measurements, since the 
redshift distribution of the AGN tends to peak at higher z than 
that for the galaxies. To account for this, for each AGN sam- 
ple we assign weights to the AGES galaxies and the random 
galaxy sample, so that their (weighted) redshift distribution is 
the same as that of the AGNs in bins of size Az = 0. 1 . Includ- 
ing these redshift weights, we rederive the autocorrelation of 
the AGES main AGES sample, correspondingly weighting the 
random galaxy sample, and multiplying each galaxy-galaxy, 
galaxy-random, or random-random pair by these additional 
weights. We then use this redshift-matched autocorrelation 
function when calculating the bias of AGNs relative to galax- 
ies. 

Furthermore, for AGNs with extended optical counterparts, 
it is interesting to compare the AGN-galaxy cross-correlation 
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with that measured for galaxies with similar color classifica- 
tion, redshift, and absolute magnitude. To create this com- 
parison sample, we assign additional weights to the AGES 
main sample galaxies, so that their (weighted) distribution in 
color, absolute magnitude, and redshift is the same as that for 
the AGN hosts (in bins of 0.25 mag, 0.25 mag, and 0.1 re- 
spectively). We then calculate the cross-correlation Wp(rp) be- 
tween this comparison sample and all AGES galaxies (again 
multiplying each galaxy-galaxy or galaxy-random pair by 
these additional weight) and derive the bias (on scales 1- 
10 and 0.3-1 /i"' Mpc) for the AGNs relative to this cross- 
correlation function. 

6.3. Absolute bias and dark matter halo mass 

As discussed in Section[T] the mass of the parent dark mat- 
ter halo may be a key parameter in determining the evolu- 
tion of galaxies and their central black holes. We therefore 
use our clustering results to estimate Mhaio for our different 
subsets of galaxies and AGNs. The masses of dark matter 
halos are reflected in their absolute bias relative to the dark 
matter distribution. To determine absolute bias (following 
[Coil et al. 2008) we first calculate the two-point autocorrela- 
ti on of dark matter a s a function of redshift, using the code 
of lSmithet al](l2003b . and assuming our standard cosmology, 
with cTg = 0.9 and the slope of the initial fluctuation power 
spectrum, F = Vl„,h = 0.21. We derive Wp{rp) as a function 
of redshift for the dark matter, integrating to tt = 25 Mpc 
as for the data. We assume an uncertainty in the model dark 
matter correlation function of 5% at each scale. 

For each subset of AGNs, we next divide the projected au- 
tocorrelation function for galaxies that are matched in redshift 
(Section l6^ by the Wp(rp) for the dark matter, averaged over 
the redshift distribution of the sample. The mean ratio be- 
tween these two correlation functions on scales 1 < < 10 
Mpc gives bg^i, where bg^^i is the absolute linear bias of 
AGES galaxies with given redshift distribution. The absolute 
bias of the AGN sample is then /jabs = ^gai^reh where bi-^i is 
the bias relative to AGES galaxies, determined by the cross- 
correlation (Section |6^ . 

Finally, we use Zjabs to derive the characteristic mass of the 
dark matter halos hosting each subset of galaxies or AGNs. 
The halo masses presented in this paper represent the virial 
mas s in the sense o f a top-hat spherical collapse model (see 
e.g..' Peebleslll993l and refer ences therein). To obtain halo 
mass, we first convert /jabs to the quantity = 5c/criM), where 
a{M) is the linear theory rms mass fluctuation in spheres of 
radius r = (3M/47rp)'/-' (p is the background density) and 
(5c ~ 1.69 is the critical overdensity r equired for collaps e. We 
convert /jabs to v using Eqn. (8) of ISheth et al.' (2001), and 
use V t o derive Mhaio following Appendix A of van den Bosch 
2002h. If we use a different relation between /^abs and v 



1 and 10 Mpc) to the Wp{rp) from the dark matter model 
of ISmith et alJ (l2003[) [calculated using the same cosmology 
assumed bv iPadmanabhan et alj (|2008|) 1. For a given 1-10 
/i"' Mpc bias, the corresp onding Mhaio from the simulations 
matches the prediction of ISheth et alJ 12001) to better than 
0.08 dex. This confirms that measuring bias from 1-10 
Mpc scales allows sufficiently accurate estimates of Mhaio- Fi- 
nally, we check that our estimates of Mhaio are relatively in- 
sensitive to our choice of erg. Decreasing trg results in weaker 
dark matter clustering and so larger ^Jabs for galaxies with a 
given observed Wp(rp). However, decreasing erg also reduces 
the Mhaio corresponding to a given b-^^s, so there is little change 
in our estimates of Mhaio- If we change cg from 0.9 to 0.8, our 
Mhaio estimates for AGES galaxies and AGNs vary by less 
than ±0.1 dex. 

In addition to top-hat virial mass, another commonly used 
mass estimator is M200, which is equal to the mass contained 
within the radius where the overdensity is 200 times the crit- 
ical density. For our chosen cosmology, M 200 is slig htly 
smaller than the top-hat virial mass, by w20% (IWhitdHoOll) . 



6.4. Uncertainties 

We determine uncertainties in Wp(rp) directly from the data, 
in two ways. For the galaxy autocorrelations, errors are de- 
rived from jackknife resampling, which generally approxi- 
mates the var iance in Wp(r „) derived from simulated galaxy 
catalogs (e.g.. lZehavi e"taLir2005l:ICoil et aLH2008l) . We divide 
the sample into five separate regions (large enough to sample 
all the appropriate scales) and recalculate Wp(rp) excluding 
each region in turn. The uncertainty in Wp(rp) is taken to be 
the variance between the jackknife samples. However, when 
calculating Wp{rp) for the AGN samples, jackknife resampling 
involves only a small number of sources, so the uncertainties 
on different scales vary significantly. Therefore, for the opti- 
cally extended AGNs we instead create a set of 100 random 
comparison galaxy samples matched in col or, a bsolute mag- 
nitude, and redshift (as described in Section l672b . but contain- 
ing only one galaxy corresponding to each AGN. We calculate 
Wp{rp) for each of the 100 samples, and take the variance to 
be the uncertainty in Wp{rp). These uncertainties are typically 
~2 times the corresponding jackknife errors. We similarly 
calculate the variance in the bias (relative to all galaxies, and 
to the matched sample) and the best-fit power law parameters, 
to account for covariance in different bins of r^. For opti- 
cally unresolved AGNs, we cannot create a matched sample 
of galaxies, so we simply repeat the calculation for the op- 
tically extended sources of each type (radio. X-ray, and IR), 
but vary the sample size to match the corresponding number 
of unresolved AGNs and take as the uncertainties the variance 
in Wp and bias between the 100 samples. 
We perform power-law fits to Wp (using the uncer- 

tainties derived from the above resampling) over the range 

Tinker et"aLl l2005h . we obtain estimates for Mhaio that are 0.3 < rp < 10 Mpc, from which we obtain the best-fit ro 

and 7. We estimate the uncertainties in ro and 7 by the vari- 
ance in the best-fit parameters between the jackknife samples 
(for the galaxy autocorrelation) or the 100 comparison sam- 
ples (for the AGN-galaxy cross-correlation), to account for 
covariance between bins of rp. We note that our jackknife 
resampling appears to underestimate the uncertainties in the 
fit parameters by about 3% of the values of ro and 7, com- 
pared to variance in simulated galaxy catalog s with compara - 
ble source densities and number of objects ( Coil et al1 l2008'). 
We therefore include a systematic uncertainty of ±3% on the 
best-fit parameters and on the absolute bias. This gives total 



similar, although slightly larger by 0.1-0.2 dex. 

We note that to estimate Mhaio, we have averaged the abso- 
lute bias on scales 1-10 Mpc. On these scales, the dark 
matter and galaxy correlation functions can have somewhat 
differ ent shapes, so that the relati ve bias is a function of scale 
(e.g., ICresswell & Percivalll2009h . To verify our method for 
estimating Mhaio using th e 1-10 /;"' Mpc bias, we compared 
the bahs-Mhaio relation of ISheth et a l. (2001 ) to a relation de- 
rived from the dark matter simulations of IPadmanabhan et al.l 
(I2b08i) . Using the simulation results for z « 0.5, we divided 
the Wp(rp) for halos of different masses (evaluated between 
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Fig. 10. — Top panel shows galaxy autocoirelations in the AGES galaxy 
sample. The points show the projected autocon'elation function in the redshift 
interval 0.25 < z< 0.8 for all galaxies in the main AGES sample (black), red 
AGES galaxies (red), and blue AGES galaxies (blue). Uncertainties are taken 
from jackknife resampling, and lines show power-law fits to the data. The 
bottom panel shows bias relative to the main galaxy sample, defined as the 
square root of the ratio between Wp(rp) for the autocorrelation of red (or blue) 
galaxies and Wp(rp) for the autocorrelation of all AGES galaxies. Real-space 
fit parameters and average bias values are given in Tabled 
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Fig. 11. — Projected cross-correlations between AGNs and all AGES 
galaxies. The top panel shows the projected cross-con'elation between ra- 
dio AGNs (orange circles). X-ray AGNs (green stars), and IR AGNs (red 
squares), and all AGES galaxies in the redshift interval 0.25 < z < 0.8. For 
comparison, the dotted lines show galaxy autocorrelations for the main AGES 
galaxy samples, weighted to have similar weighted redshift distributions to 
the combined sample of AGNs. The bottom panel shows bias relative to 
AGES galaxies, defined as the ratio between Wp(rp) for the AGN-galaxy 
cross-correlation and Wp(rp) for the galaxy autocorrelation. 



uncertainties for the galaxy autocorrelations that are similar 
to those derived from simulations, but has a negligible effect 
on the (much larger) errors for the cross-correlation of AGNs 
and galaxies. 

7. CORRELATION RESULTS 

In this section we discuss the results of the correlation anal- 
ysis and the characteristic dark matter halo masses for galax- 
ies and AGNs. 

7.1. Galaxy autocorrelation 

As a first exercise, we measure the autocorrelation of the 
main AGES galaxy sample, as well as for the red and blue 
galaxy samples separately. The resulting Wp(rp) for the galaxy 
autocorrelations are shown in Figure [TO] and the parameters 
of power-law fits to the data (including the reduced x^) are 
given in Table |2l The range of absolute magnitudes for the 
red and blue galaxy samples are -19 < MoA,.-5logh < -23 
and -19 < Moi,.-51og/i < -22.5, respectively. The autocor- 
relation of the main galaxy sample is best fit by ro = 4.6 ± 0.2 
/i"' Mpc and 7 = 1.8 ± 0.1. This is bro adly consistent with 
previous re sults [see, e.g.. Table 4 in Brown et aTl (l2003h . 
Table 1 of ICoil et all (|2008|) . iBrown et al. (200^and ref- 
erences therein]. The corresponding absolute bias indicates 



that typical AGES galaxies reside in dark matter halos with 

Dividing the galaxy sample into red and blue subsets, we 
find that red galaxies show a significantly stronger autocor- 
relation and a steeper slope than the blue galaxies, reflecting 
the well-established relationship between color and cluster- 
ing. The difference in clustering between red and blue galax- 
ies is not simply a luminosity effect. If we select random 
samples of red and blue galaxies that are uniformly distributed 
in bins of 0.25 mag between -19 > Mo i^ > -21.25, the rela- 
tive bias between red and blue galaxies changes by only 2%. 
We note that the Wp(rp) is not exactly a power law, especially 
for the red galaxy sample. The apparent upturn of Wp(rp) on 
scales less than 1 /;"' Mpc is indicative of an additional "one- 
halo" term in the correlation function, corresponding to dark 
matter halos that contain multiple galaxie s (e.g., Zehavi et aD 
'2004; Zheng et alj|2003 ICoil et alJlIOOSt IBrown et al.ll2008r 
iZheng et al...2009l) . 

As a check, we also perform the galaxy autocorrelation 
using an alternative clustering measure, the oj statistic of 
|Padmanabhan et al. (2007). This statistic has the advantage of 
being less sensitive than Wp(rp) to large-scale structures and 
the effects of truncating the line-of-sight integral at Tr^ax- In 
addition, it is possible to evaluate uj{R) by converting the inte- 
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TABLE 2 
Correlation results^ 
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0.3 


1.05 ±0.12 0.65 ±0.22 


1.01±0.11 


0.66 ±0.23 


1.34±0.15 


12.8:1 


IRAC 


133 


0.42 


3.1 ±0.7 


1.2±0.3 


0.3 


0.73±0.18 0.26±0.33 


0.70±0.17 


0.34 ±0.36 


0.90 ±0.20 
















AGNs (unresolved)^ 








X-ray 


121 


0.57 


5.0±0.7 


1.7 ±0.2 


1.3 


1.14±0.17 0.69±0.25 






1.52±0.23 


13 0+"^ 


IRAC 


105 


0.55 


4.1 ±0.7 


2.0 ±0.2 


2.2 


0.70±0.18 0.81 ±0.27 






0.93 ±0.24 


11.61 



Note. — All units for rp and ro are in li ' Mpc (comoving). 

^ Results in the top three rows are for the autocoirelation of all AGES main sample galaxies, red AGES galaxies, and blue AGES galaxies, derived using the lLandv & Szala^ 
119931) estimator (Eqn.l4l. Values in the other rows are for the cro.¥.9-con"elation of the given AGN sample with all AGES main sample galaxies, derived using the estimator given 
by Eqn. |6).^ Median redshift for the sources in the sample.'^ Bias relative to AGES main sample galaxies with matched redshift distributions (§ 16. 2L averaged over the scales 
shown.'' Bias relative to AGES main sample galaxies with matched distributions in color, absolute magnitude, and redshift (§ 16. 21 , averaged over the scales shown. ^ Virial mass in 
the sense of a top-hat collapse model, assuming a cosmology with n,„ =0.3, Oa = 0.7, and crs =0.9.^ AGNs selected both in the X-rays and IR.^ Optically unresolved radio AGNs 
are not included in the clustering analysis, owing to the small sample size and corresponding large uncertainties. 



grals to Riemann sums, and so it does no t require binning the 
data. Using the (2,2) filter preferred in iPadmanabhan et al.l 
(120071) . and evaluating their Eqn. (29), we derive lo(R) and 
calculate jackknife uncertainties as described above. We per- 
fo rm a power-law fit to uj and derive ro and 7 using Eqn. (24) 
of lPadmanabhan et al.l (l2007h . The best-fit values are ro = 4.8 
and 7 = 2.0, within 5% and 7%, respectively (^Ict) of the 
values derived from the Wp{rp) fit. We note, that since u:{R) 
involves an integral over rp as well as tt, the bins of ojifi) 
are not independent. We therefore use Wp{rp) to estimate the 
AGN-galaxy cross-correlation and relative bias as a function 
of scale. 

7.2. AGN-galaxy cross-correlation 

We next measure the cross-correlation of AGNs with AGES 
galaxies, which gives us the bias of AGNs relative to galax- 
ies, from which we derive the absolute bias and Mhaio for the 
AGNs. We first perform the AGN-galaxy cross-correlations 
for all AGNs (including those with extended or unresolved 
counterparts). The Wp(rp) values are shown in Figure[TT] and 
the parameters of power-law fits to the correlation functions 
are given in Table |2l We further divide the AGN samples 
into those with extended or unresolved optical counterparts 
(Figs. [T2I and [T3]l. For optically extended sources, we com- 
pare the clustering of AGN hosts to a control sample of nor- 
mal galaxies that are appropriately weighted to have similar 
colors, absolute magnitudes, and redshifts as the AGN hosts 
(Section I6.2l i. Here we discuss the results for the different 
classes of AGN, which show significant differences in clus- 
tering. 

7.2.1. Radio AGNs 

The radio AGNs are most strongly clustered, with signifi- 
cant bias relative to normal galaxies indicating that they re- 
side in massive dark matter halos (Mhaio ~ 3 x 10'^ /i"' Mq) 
characteristic of large galaxy groups or small clusters. There 



is an increase in bias on scales less than 0.5 h Mpc, simi- 
lar to that seen for red AGES galaxies. This suggests that on 
small scales, the cross-correlation of radio AGNs with galax- 
ies is dominated by pairs of objects in the same dark matter 
halo, which is to be expected if these sour ces reside in massive 
systerns (e.g.,|Ze havi et al. 2004; Zheng e t alJ2007l:[Coil et alJ 
l2008HBrowneraL.2008: .Zheng et aL.2009h . 

Interestingly, compared to a control sample of AGES galax- 
ies matched in color. Mo 1^, and z, the radio AGNs show no sig- 
nificant difference in clustering on all scales from 0.3 to 10 
Mpc. This is apparently in contrast to recent studies that found 
that radio AGNs are more strongly clustered (and have corre- 
spondingly larger halo masses) co mpared to matched contro l 
samples of galaxies. For example, Mandelbaum et al.] (l2009l) 
selected radio AGNs from SDSS and found that their dark 
matter halos are roughly twice as massive as those for galax- 
ies matc hed in redsh i ft, stel lar mass, and stellar velocity dis- 
persion. IWake et al] (l2008l) found a similar effect for more 
luminous radio AGNs at 0.4 < z < 0.7 taken from the 2SLAQ 
Luminous Red Galaxy survey, with the radio sample inhabit- 
ing ^ 2 times larger halos than the control galaxies matched 
in luminosity and color. 

We note that the differences between these studies and our 
own results are only marginally significant; if the bias of the 
Bootes radio AGNs were increased by 1 .5(7, this would corre- 
spond to a factor of 2 larger halo mass re lative to the control 
gala xies, and w ould match the results of iMandelbaum et alJ 
and IWake et alJ N onetheless, it may be worth noting that 
IMandelbaum et all studied a subsample of high stellar-mass 
objects that corresponds closely to the range of Pi.4ghz and 
host galaxy mass probed in Bootes (assuming mass-to-light 

ratios from iKauffmann et al 2003b,) . and found that for these 

objects, the difference in clustering between the radio AGNs 
and the control sample is small. 

These results may hint at an interesting point regarding the 
fueling of radio AGNs. In particular, the importance of envi- 
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Fig. 12. — Same as Figure fTTI for only those sources with extended optical 
counterparts (Section [5). Real-space fit parameters and average bias values 
are given in Table[2] 



ronment in triggering radio activity may depend on the radio 
luminosity relative to the galaxy's stellar mass. Compared 
to the Bootes radio AGNs, the full Mand elbaum et al.. sample 
extends to lower stellar masses, while the IWake et al l AGNs 
have similar host galaxy absolute magnitudes but higher ra- 
dio luminosities (the ir median fi4GHz is ^^3 times higher 
than in our sample; ISadler et"al1 12007|) . Therefore, both 
samples probe sources with a higher typical ratio of radio 
power to stellar m ass than the Bootes sample (or in the 
iMandelbaum et all high-mass sample), for which we see little 
clustering difference between radio AGNs and control galax- 
ies. 

This suggests that for galaxies with low radio luminosities 
(relative to their stellar masses), environment may have little 
effect on the fueling of radio activity. This is certainly true in 
the limit of very low luminosity rad io AGNs, which ar e found 
in virtually every massive galaxy (Sadl er et al.lll989l) . How- 
ever, the fueling of AGNs with higher radio power relative 
to galax y mass, which ar e less common than the low-power 
sources ( lBestetalJl2005h . may be more strongly dependent 
on the mass of the host dark matter halo. In the future, a 
larger study that examines radio AGN clustering over a wide 
range of radio luminosities and environments will allow us to 
better understand how radio activity is triggered in halos of 
different masses. 

7.2.2. X-ray AGNs 

For the X-ray AGNs, the clustering is consistent with that of 
typical AGES galaxies on scales of 1-10 Mpc, indicating 
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Fig. 13. — Same as Figure fTTI for only those sources with unresolved op- 
tical counterparts (Section [3). There are only nine radio-selected, optically 
unresolved AGNs, so the corresponding Wp(rp) values have large uncertain- 
ties and are not included here. Real-space fit parameters and average bias 
values are given in Table[2] 



that they reside in halos with Mhaio ~ 10 Mq, character- 
istic of poor to moderate-sized galaxy groups. There is little 
difference in the clustering between the optically unresolved 
and and optically extended X-ray AGN samples. This sug- 
gests that obscured or optically faint AGNs show similar clus- 
tering to their optically bright counterparts. This is analogous 
to the well-established results for optical quasars, for which 
clustering is only weakly (if at all) dependent on optica l lumi- 
nosity (e.g., Adelberger & Steidel 2005; Porciani & N orberd 
l2006tlMvers et alj|2007l;[da Angela et al.ll2008l) . although un- 
like in optical quasar samples, some of our X-ray AGNs are 
likely optically faint because of obscuration, rather than low 
intrinsic luminosity. Relative to a matched sample of normal 
galaxies, the optically extended X-ray AGNs show no signifi- 
cant bias. 

However, on scales less than 1 /i"' Mpc, the clustering of 
X-ray AGNs differs significantly from that of normal galax- 
ies. The slope of the AGN-galaxy cross-correlation function 
for X-ray AGNs is significantly flatter than the autocorrela- 
tion function for normal AGES galaxies, such that the X-ray 
AGNs are significantly antibiased on small scales (0.3-1 
Mpc). This antibias on small scales is found for all X-ray 
AGNs relative to all AGES galaxies, and also for optically ex- 
tended X-ray AGNs, relative to the control sample of normal 
galaxies with matched optical properties. 

A similar result was recently obtained by iLi et alj (l2006h 
for optically selected narrow-line AGN in SDSS. Specifically, 
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iLi et alj (l2006l) found that on scales greater than 1 Mpc, 
the cross-correlation between optical AGNs and galaxies is 
similar to a control sample of normal galaxies, but that the 
AGN-galaxy clustering is weaker on scales less than 1 /i"' 
Mpc. Comparing to mock galaxy catalogs taken from cos- 
mological simulations, Li et al- (2006) showed that the small- 
scale antibias can naturally be explained if the AGNs are pref- 
erentially located in central galaxies within their dark matter 
halos, which would serve to suppress the "one-halo" term cor- 
responding to strong clustering on small scales. While statis- 
tical uncertainties in the clustering analysis prevent us from 
performing a similar sophisticated comparison here, the con- 
sistent antibias on small scales for the X-ray AGNs suggests 
that they may also preferentially reside in central galaxies. 

7.2.3. IRAGNs 

For the IR-selected AGNs, the cross-correlation with AGES 
galaxies is weaker and has a flatter slope than for radio or 
X-ray AGNs. The estimated dark matter halo masses are 
Mhnio < 10'^ Mq. These halo mass estimates are un- 
certain because at these redshifts and halo masses, b^bs is a 
weak function of Mhaio- However, we can conclude that IR 
AGNs typically reside in halos that are significantly smaller 
than those that host radio or X-ray AGNs. We also note that 
the reduced for the fit to Wp(rp) is quite low (0.4) for the IR 
AGNs, possibly because of covariance between the Wp(rp) on 
different scales that is not reflected in the errors on the indi- 
vidual points. The uncertainties in the fit parameters and the 
bias are derived directly from resampling, and so account for 
any covariance between scales. 

In contrast to the radio or X-ray AGNs, optically extended 
IR AGNs are significantly antibiased relative to normal galax- 
ies with similar optical properties. This suggests that the pro- 
cess that triggers IR-bright AGN activity depends not only 
on host galaxy properties but also on local environment, with 
AGNs preferring underdense regions. The optically extended 
IR AGNs are less biased (more strongly antibiased) on small 
scales (0.3-1 /i"' Mpc) compared to larger scales (1-10 /i"' 
Mpc), reflected in the flat slope of the cross-correlation func- 
tion. Keeping in mind that the relative bias on the two scales 
differ by only 1.5cr, this may suggest that IR AGNs preferen- 
tially reside in the central galaxies of their (relatively small) 
dark matter halos. 

7.2.4. X-ray and IR AGNs 

Finally, we examine the AGN-galaxy cross-correlation for 
objects that are selected as both X-ray and IR AGNs. Owing 
to the limited sample size, we do not divide these sources into 
optically extended and unresolved subsets. The clustering pa- 
rameters, bias, and characteristic Mhaio for these sources are 
given in Table|2](to avoid confusion, the cross-correlation re- 
sults for this sample are not shown in Figure [TTli. AGNs that 
are selected in both the IR and X-rays show essentially iden- 
tical clustering to the full sample of IR-selected AGNs, with 
significant antibias relative to AGES galaxies. This suggests 
that the IR-selected AGNs are weakly clustered, regardless of 
whether they are selected in X-rays, and are found in smaller 
dark matter halos than typical X-ray AGNs (most of which 
are not selected with IRAC). We discuss the implications of 
the clustering results for all the AGN samples in Section|9] 

8. X-RAY SPECTRA AND EDDINGTON RATIOS 

In this section, we probe the properties of the accre- 
tion process in the different classes of AGNs, by measur- 



ing average X-ray spectra and Eddington ratios. There is 
growing evidence that black hole accretion can occur in 
different states, with a key paramete r being the Edding- 
ton ratio (e.g., [Ch urazov et al. 2005; Me rloni & Heinzl 120081: 
Kauffmann & Heckm an 2009), and AGNs in these different 
modes may play distinct roles in galaxy evolution. In addi- 
tion, many AGNs are obscured by gas and dust; according 
to "unified models", this obscuration originates in a roughly 
toroidal structure that is intrinsic to the accretion flow (e.g., 
lAntonuccill 19931: lUrrv & Padovanilll995b . We first compute 
average X-ray spectra, which can be used to detect the ab- 
sorption by gas or can possibly indicate different accretion 
modes. We then derive bolometric AGN luminosities (Lboi) 
for the different classes of AGNs, which we use along with 
the estimates of black hole masses to derive Eddington ratios. 

8.1. X-ray stacking analysis 

To estimate the average X-ray spectra for AGNs, we use 
an X-ray stacking analysis. Unabsorbed, Type 1 AGNs typi- 
caUy have power-law X-ray spectra with F ~ 1 .8 (iTozzi et al.l 
l2006l) . but intervening neutral gas preferentially absorbs soft 
X-rays and so hardens the X-ray spectrum. Alternatively, 
a harder X-ray spectrum may indicate a different accre- 
tion mode than is typical for optically bright Seyferts or 
quasars. Galactic black hole binary systems show multi- 
ple X-ray spectral states with varying spectral hardness that 
are caused by changes in the accr e tion mode (for reviews 
see 'R emiflard & McClintocfl [2006'; Naraya n & McClintocfl 
,2008i) . The "low-hard" state found in black hole binaries is 
often interpreted as a radia tively inefficient mod e (lEsin et al.l 
Il997h . such as an ADAF (iNarayan & Yilll995h . Unusually 
hard X-ray spectra observed in an AGN (particular for a low 
accretion rate) may indicate a similar mode of accretion (e.g., 
Hopkins et al. 2009b), as has been sugges ted for very low E d- 
dington AGNs in red galaxies in AGES ( Brand et al.ll2005l) . 

Even for the X-ray-detected sources in our AGN sample, 
we cannot accurately determine spectral shapes for individ- 
ual sources since most have <10 counts. However, we can 
use X-ray stacking to calculate average spectral properties for 
subsets of sources. Specifically, we calculate the hardness ra- 
tio, defined as 

where H and S are the observed counts in the hard (2-7 keV) 
and soft (0.5-2 keV) bands, respectively. 

The stacking pro cedure is described in Section 5.1 of 
iHickox et al.l (l2007h . Stacked counts are defined as the num- 
ber of (background-subtracted) photons detected within rgo 
from the source position, where rgo is an approximation of 
the 90% point-spread function (PSF) energy encircled radius 
at 1.5 keV, and varies as": 

r9o= 1"+1O"(0/1O')^ (13) 

We subtract a small background of 3.0 count s~' deg"^ 
source"' in the 0.5-2 keV band and 5.0 count s"' deg"^ 
source"' in the 2-7 keV band (note that we do not include 
sources found in 1 1 pointings with high flaring background; 
see Section 5.1 of Hickox et al. 2007). In order to maximize 
the number of source counts, we do not limit the stacking to 
the central 6' around the pointing center for each observation 

'* Chandra Propose r's Observatory Guide 

|http://cxc.harvar d.edu/pro poser/POG/| 
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TABLE 3 

X-RAY HARDNESS RATIOS AND SPECTRAL SHAPES 







A 11 

All 


< 20 counts 


Subset 




HR 


HR 








All AGNs 






Radio 


103 


-0.44 ±0.04 1.6 ±0.1 


-0.42 ±0.10 


1.5±0.2 


X-ray 


318 


-0.34±0.01 1.4±0.1 


-0.24 ±0.02 


1.1±0.1 


IR 


198 


-0.34 ±0.02 1.4 ±0.1 


-0.11 ±0.04 


0.9±0.1 






AGNs (galaxies) 






Radio 


95 


-0.37 ±0.06 1.4 ±0.1 


-0.36 ±0.11 


1.4±0.2 


X-ray 


212 


-0.21 ±0.02 1.1 ±0.1 


-0.15 ±0.03 


1.0±0.1 


IR 


114 


-0.10 ±0.04 0.8 ±0.1 


0.11 ±0.06 


0.4 ±0.1 






AGNs (unresolved) 






Radio 


8 


-0.51 ±0.06 1.7 ±0.1 


-0.71 ±0.27 


2.2 ±0.6 


X-ray 


106 


-0.46 ±0.02 1.6 ±0.1 


-0.41 ±0.03 


1.5±0.1 


IR 


84 


-0.44 ±0.02 1.6 ±0.1 


-0.31 ±0.05 


1.3±0.1 



^ Number of sources in stacked sample; excludes sources in Chandra exposures with high 
background.^ Effective power law photon index, derived from HR using Eqn. ll4l 



(unlike in lHickox et alj|2007h . Count uncertainties are calcu- 
lated using the approximation <Jx = VX + 0.75+1, wh ere X 
is the number of counts in a given band (lGehrelslll986l) . Un- 
certainties in HR are derived by propagating these count rate 
errors. 

To estimate an effective power-law photon index F, we use 
PIMMS to calculate the HR that would be observed for dif- 
ferent values of F, given the ACIS Cycle 4 on-axis response 
function. For the small Galactic A^h toward this field (~10^° 
cm~^), the relationship between HR and F can be approxi- 
mated by 

HR = 0.29-0.46F. (14) 

It is possible that a few very bright sources can dominate 
our estimates of the average HR, so we also calculate HR 
for only those sources with less than 20 counts. The aver- 
age HR values and the corresponding F estimates are given 
in Table |3] As a check, we perform the same analysis for 
the X-ray-detected sources using the observed counts in the 
iKenter et al.l (l2005[) catalog [r ather than the full stacking pro- 
cedure of iHickox et al.l (|2007|) 1. and obtain essentially identi- 
cal results. 

As shown in Table [3] X-ray and IR AGNs with unresolved 
optical counterparts have relatively soft average spectra, with 
F ~ 1.6. The bright optical emission in these sources sug- 
gests that the nucleus is relatively unobscured by dust, and 
their X-ray spectra are typical of X-ray AGNs that have little 
intervening absorbing gas. Optical obscuration and X-ray ab- 
sorption are strongly correlated for most AGNs (iTozzi et al.l 
I2OO6) . so these sources are consistent with being a popula- 
tion of unobscured Seyfert galaxies. In contrast, the X-ray 
and IR AGNs with extended optical counterparts have sig- 
nificantly harder average X-ray spectra, with 0.8 < F < 1.4. 
These hard X-ray spectr a may be caus ed by absorption by gas, 
with A^H ~ 10^^ cm"2 (iHickox et all 12007) . or alternatively, 
these sources might have intrinsically hard X-ray spectra ow- 
ing to radiatively inefficient accretion. In either case, it would 
be expected that the sources would have relatively weak op- 
tical emission, as observed. We conclude that the average X- 
ray spectra of the X-ray and IR AGNs are consistent with the 
properties of their optical counterparts. 

Among the radio-selected AGNs, the six sources also de- 
tected in X-rays have relatively soft spectra with HR ~ -0.4, 
corresponding to a typical AGN spectrum with little X-ray ab- 
sorption. However, the radio AGNs that are not individually 
detected in X-rays have a softer stacked spectrum; the aver- 



age number of background-subtracted counts per source are 
0.41 ± 0.08 in the 0.5-2 keV band and 0.10 ± 0.04 in the 2-7 
keV band. If all this flux comes from the AGN, it corresponds 
to a soft spectrum with F ~ 1 .9. The lack of hard X-ray flux 
in these sources provides an upper limit to the contribution 
from an AGN, assuming F = 1.8. The average AGN flux for 
these radio galaxies must be 5o.5-7keV < 2.1 x 10"'^ erg cm"^ 
s"' (3 a upper limit). For the range in redshifts observed for 
the radio AGNs, this corresponds to an upper limit on Lx of 
(0.5-5) X 10^2 ergs-i. 

8.2. Eddington ratios 

We further use the observed luminosities to estimate Ed- 
dington ratios, which can provide insights as to the accretion 
states of the AGNs. For a central object with mass Mbh, 
the E ddington luminosity is given by ( Shapiro & Teukolskvl 
IT983h 

(15) 



LEdd = 1.3 X 1O^''*(Mbh/M0) erg s" 



The Eddington ratio is defined as A = Lboi/i'Edd, where Lboi 
is the bolometric accretion luminosity of the system. Lu- 
minous quasars tend to have relatively high A > 0.1 (e.g, 
McLure & Dunlop 2004; Kollmeier et al. 2006), but some 
Seyfert galaxies, LINERs, and low-luminosity X-ray and ra- 
dio AGNs can have A < 10"^, and even as low as A ~ 10"^ 
(e.g.. lHoll2002l:ISoria et aUlIOOaiSikora et al.ll2007^ . 

We calculate a rough estimate of A for the different classes 
of AGNs with extended optical counterparts (for which we 
can estimate the galaxy stellar mass, and thus Mbh)- To es- 
timate Lboi for X-ray AGNs, we scale from the observed Lx- 
For radio AGNs that are not detected in X-rays, we use the 
X-ray stacking results (Section lSTt to derive approximate up- 
per limits on the AGN luminosity. We take the 3(7 upper limit 
on the average nuclear 0.5-7 keV flux, and derive the corre- 
sponding upper limit on the luminosity for each source. 

We next obtain Lboi by scaling from Lx- We first convert 
Lx from the 0.5-7 band to the 0.5-8 keV band by multi- 
plying by 1.1, which is roughly valid for power-law spec- 
tra with photon index 0.8 < F < 1.8. We then convert to 
Lboi using the lumino sity-dependent bolometric corrections of 
iHopkinsetall (l2007h . which are in the range BC(),5-8 keV — 
10-20. This correction assumes that the intrinsic broadband 
SED of these AGNs is similar to that of unobscured quasars. 
If instead the source has an SED that is more dominated by 
X-ray emission, this corre ction may overestimate Lboi (e-g-, 
IVasudevan & Fabianll2007h . The distributions in Lboi for the 
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Fig. 14. — Distributions in (a) L|,oi, (b) Mbh, and (c) Eddington ratios for 
AGNs with extended optical counterparts, for radio. X-ray, and IR AGNs. 
The Li,oi and Lbol/^Edd estimates for the X-ray undetec ted r adio AGNs show 
onl y upp er limits derived from X-ray stacking (Section fS.lt . See text in Sec- 
tion [Ojfor details. Radio, X-ray, and IR AGNs have progressively smaller 
typical black hole masses and higher Eddington ratios. 



X-ray sources (and limits for the radio AGNs not detected in 
X-rays) are shown in Figure [T4b. 

For IR AGNs that are not detected in X-rays, Lboi is de- 
rived from th e rest- frame 4.5 /im luminosity. Using the 
iHopkins et al.l (120071) scaling somewhat overestimates Lboh 
compared to the values derived from the X-rays (for those 
sources that are both X-ray and IR AGNs). This may be ow- 
ing to the fact that the host galaxy can contribute some of the 
observed IR flux. We therefore decrease the IR bolometric 
corrections by a factor of 2 so that Lboi derived from the X- 
rays and IR are consistent, and use these corrections for all 
the IR AGNs. The distributions in Lboi for the IR AGNs are 
also shown in Figure [T4b. 

For an estimate of Mbh, we first determine the bulge lumi- 
nosity in the B band (Lg bui) a nd then scale, using the rela- 
tion of iMarconi & Hunt! (l2003h . between 

i'B.bui and Mbh ■ We 

derive the absolute B band magnitude of the galaxy, Mb, by 
adding the rest-frame ° ' (§~ '") color to the evolution-corrected 
absolute magnitude in the *' 'r band (Section l3.3l l, and then 
subtracting 0. 1 mag for the small conversion between the '^■^g 
and B bands. We then calculate Lb (in Lq), taking the absolute 
magnitude of the Sun in the B band to be 5.48. 

We further require an estimate of the ratio of bulge to to- 
tal luminosity in the B band, {B/T)b, which can vary from 
less than 0. 1 for Scd galaxies to 1 for ellipticals. Here we use 




log(M3H [Md) 

Fig. 15. — Black hole mass funct ion f or 0.25 < < 0.8, derived from Mbh 
measurements described in Section lS^l compared to local mass function es- 
timates in the literature. The thick black line shows </<(Mbh) derived from the 
AGES data, including rough completeness corrections. The sohd cyan and 
d ashed green hues show local estimates of (/>(Mbh) compiled in Figure 11 
of lLauer et alj 12007), derived from luminosity functions and distributions in 
bulge stellar velocity dispersions, respectively. Our rough estimates of A/bh 
give a BH mass function that is similar to those derived from other methods. 



the morphologi cal analysis of the Millennium Galaxy Cata- 
log (MGC) by Allen et all (I20b6). who calculated structural 
parameters (including bulge and disk luminosities) for a large 
sample of galaxies. Using a subset of 8243 galaxies from the 
MGC with -24 < Mb < -18 (roughly the magnitude disti-i- 
bution of our AGN host galaxies), we calculate the average 
(B/T)b in bins of M-r color. {B/T)b varies from ~0.1 at 
M-r = 1 (the bottom of the blue cloud) to ~ 0.5 at M-r = 3 (the 
top of the red sequence), with variance of A(B/T)b ^ 0.2 in 
each color bin. Using this relation and the rest-frame ^-^(u-r) 
colors (making a slight correction of 0.1 mag to convert to 
u-r), we derive (B/T)b, and thus Ls.bub for the AGN hosts. 

Finally, we derive Mbh from Lg bui using the relation of 
iMarconi & Hunll (1200 3) for their "Group 1" sample of galax- 
ies. The corresponding distributions in Mbh and A are shown 
in Figs. \T4h and[T4fc. Owing to the scatter in the above re- 
lations, these Mbh estimates may be uncertain by as much as 
0.5 dex for individual galaxies, but we expect the average dis- 
tribution to be accurate to within ^ 0.2 dex. The method for 
measuring Mbh likely introduces a further systematic factor 
of ^ 0.3-0.5 dex, but this does not affect the relative distri- 
butions in Mbh and A for the different samples, which are 
primarily of interest here. 

As a check on these Mbh estimates, we derive a rough 
SMBH mass function from the AGES galaxies and compare 
with those derived by previous works. For galaxies with 
Mbh ^ 1O^M0 the optical luminosity is low enough that they 
are not detected out to z = 0.8. To account for this, we use 
the observed distribution of black hole masses with redshift 
to estimate the typical volume probed for sources as a func- 
tion of Mbh. Including the galaxy sampling weights (Sec- 
tion [STjl, we derive the mass function (/)(Mbh). We directly 
co mpare this to loca l black hole mass functions compiled 
by 'La uer et aTl (l2007h . and find that they closely agree (Fig- 
ure[T5]l, suggesting that on average our BH mass estimates are 
reasonably accurate. We note that we can directly compare 
our BH mass function for objects at 0.25 < z < 0.8 to 0(Mbh) 
derived locally, because the high-mass end of the BH mass 
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function has evolved only slightly since the epochs probed by 
the t ypical AGES galaxies. The observed AG N "downsizing" 
(e.g.j Hasinger et al.l2005tlBarger et al.l2005b implies that BH 
growth at z < 0.5 is dominated by relatively small black holes, 
so that the change in the BH mass function from z ^ 0.5 to 
z = is only 0. 2 dex at M rh 10^ Mq, and negligible for 
Mbh > 10^ Mq dMerloni & Heinz 2008). 

Figure [T4jc) shows that the different classes of AGNs have 
significantly different characteristic Eddington ratios. X-ray 
AGNs have a wide range of 10"^ < A < 1, while the IR AGNs 
have higher Eddington ratios, almost all having A > 10"^. The 
X-ray undetected radio AGNs have very low upper limits on 
A, with most sources having A < 10"^. Although these es- 
timates of A are highly uncertain, particularly for individual 
galaxies, we can conclude that the population of X-ray AGNs 
extends to lower Eddington ratios than that of the IR AGNs, 
while the radio AGNs have very low Eddington ratios. 

Finally, we note that the sources that are selected as both 
X-ray and IR AGNs have a distribution in A that is almost 
identical to that of the full IR AGN population. The median 
values of A for the IR AGNs that are detected and undetected 
in X-rays differ by only 0.05 dex. We note as well that the 
clustering of the X-ray and IR AGNs is essentially identical 
to that for the full IR AGN sample (Section \1 .2A\ . and that 
the stacked X-ray flux from the X-ray undetected IR AGNs is 
consistent with having low-level X-ray emission (as shown in 
Appendix iBli. Therefore, the IR-selected sources may repre- 
sent a single coherent population of AGNs, distinct from typ- 
ical X-ray AGNs, most of which are not selected with IRAC. 

9. DISCUSSION 

In this paper we have used observations across a broad 
range of wavelengths together with spectroscopic redshifts 
to provide a general picture of the populations of accreting 
SMBHs, their host galaxies, their environments, and their ac- 
cretion modes at moderate redshifts (0.25 < z < 0.8). As dis- 
cussed in Section [T] there is evidence for a connection be- 
tween accretion onto SMBHs and the evolution of their host 
galaxies. Therefore the properties of AGN and their hosts that 
we observe at z ^ 0.5 provide a "snapshot" in the co-evolution 
of these objects over cosmological time. 

9.1. A simple model of AGN and galaxy evolution 

In this section we present a simple "cartoon" model of AGN 
and host galaxy evolution that explains, in broad terms, the 
observations described in this paper We base this model on 
four ideas that have received much discussion in the literature 
on AGN and galaxy evolution, and are supported by consid- 
erable observational and theoretical evidence: 

1 . In the early universe, galaxies initially form as systems 
that are rich in gold gas and have rotation-dominated 
dynamics (although they can be clumpier and more tur- 
bulent than present-day disk galaxies). This scenario 
is supported by observations o f high-redshift galaxies 
(e.g., Ravindranath et al. 2006t lEknepreen et al.ll2007l) 
as well as theoretical argument s for the dissipational 
collapse of cold clou ds (e.g., IWhite & ReesI 119781: 
iFall & Efstathioulll980l). and de t ailed numerical mod- 
elsje.g.. [Robertson et al.1 12004 l2006t iHopkins et all 
l2009ah . 

2. Optical quasars are found in dark matter halos with 
a characteristic mass of Mhaio ~ 3 x 10'^ Mq that re- 



mains effectively constant with redshift. In large sur- 
veys, the bias of quasars is observed to increase with 
redshift, consistent with this characteristic pare nt dark 
matte r halo mass (Porciani et al. 2004; Croor net al.l 
200l [Porc iinr&Norberg 2006; Coil et al. 200l 
Myers et al .1 12007; Shen et al. 2007: da Angela et all 
2008r iPadm anabha n et al. Quasars represent 

pha ses of accr etion at high Eddington ratios (A > 
0.1;|McLure&^Dunlop 2004; Kollmeier et al. 2006), in 
which massive black holes accrete the bulk of their final 
mass (e.g., lYu & Tremain3l200l IShankar et aT]|2009t 
lYu & Lull2008^ ■ 

3. The processes that fuel quasars are also responsi- 
ble for the creation of stellar bulges. There is ev- 
idence that the quasar phase is often preceded by 
dust-obscured phases of very rapid star formation, 
which would be manifest observationally as ultralu- 
minous infrared g alaxies (ULIRGs) or submilhmeter 
galaxies (SMGs) (S anders & M irabel'1 996l:ISmail et all 
[1997). The process that disturbs the galaxy and fuels 
the star formation and accretion can also likely cre- 
ate a dynamically hot stellar bulge. Possible physi- 
cal proces ses include major merg ers of gas-rich galax- 
ies (e.g., Kauffmann & Haehnelt 2000; Springe l et aLl 
2005 ; Hopkins et al. 2006b) and disk instabiliti es (e.g ., 
Moetalll998l: iBower et al.ll200'6l: iGenzel et al.li2008h 

4. After the creation of the stellar bulge, the star formation 
in a galaxy is quenched on a relatively short timescale. 
The quenching can occur because the available gas 
gets used up in starbursts, or because the dark matter 
halo has increased in mass so that the gas in the sys- 
tem is heated to the high virial temperature of the halo 
and cannot cool efficiently (e.g..' Rees & Ostrikerll 19771: 
iBirnboim & De kel 2003). Feedback from an AGN may 
also help to blow away gas and quench star formation 
(e.g., Hopkins et al. 2006b). Further mechanical feed- 
back (likely also from an AGN) prevents subsequent 
gas coolin g and keeps the sys tem "red and dead" (e.g., 
Churazov etal .112005 ; Bower et all 120061; ICroton et al.1 
2006tiKhalatvan et al.ii2008i) . 

These ideas imply a relatively straightforward evolutionary 
sequence for massive galaxies and their central black holes. 
Most aspects of this picture have been addressed in much 
more detail elsewhere (e.g., Kauffmann & Haehnelt' '20001 
Springel et al. 2005; Croton et al. 2006; Bower etal. 200^ 
Mona co et al.ll2007l: IHopkins et al.l i2008a b; Somerville et all 
,2008i) ; here we present only a simple overview and show that 
it is generally consistent with our results. Comparisons of the 
observations with more detailed model predictions will be left 
for future work. 

9.2. Evolutionary sequence 

In our simple picture, the sequence of AGN and galaxy evo- 
lution begins in the early niverse (z > 6), where galaxies are 
formed as gas-rich, rotation-dominated systems and contain 
small "seed" central black holes. As their parent dark mat- 
ter halos grow up to masses ^10'^ Mq, these galaxies retain 
their general morphological and dynamical properties, includ- 
ing flattened morphologies, small stellar bulges, and small 
central black holes. This growth of these galaxies occurs 
through the accretion of further cold gas or through minor 
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Fig. 16. — Schematic for a simple picture of AGN and host galaxy evolution. The picture consists of an evolutionary sequence that occurs at different redshifts 
for halo masses of different sizes. In this scenario, luminous AGN accretion occurs preferentially (thi'ough a merger or some secular process) when a host dark 
matter halo reaches a critical Mhaio between lO'^ and lO'^ /i"' Mq (this phase is indicated by the solid ovals). Once a large halo reaches this critical mass, it 
becomes visible as a ULIRG or SMG (owing to a burst of dusty star formation) or (perhaps subsequently) as a luminous, unobscured quasar The ULIRG/quasar 
phase is associated with rapid growth of the SMBH and formation of a stellar spheroid, and is followed by the rapid quenching of star formation in the galaxy. 
Subsequently, the young stellar population in the galaxy ages (producing "green" host galaxy colors), and the galaxy experiences dechning nuclear accretion 
that may be associated with an X-ray AGN. Eventually the aging of the young stars leaves a "red and dead" early-type galaxy, which experiences intermittent 
"radio-mode" AGN outbursts that heat the sun'ounding medium. For "medium" initial dark matter halos, the quasar phase and formation of the spheroid occur 
later than for the system with high halo mass, so that at z ~ 0.5 we may observe the "green" X-ray AGN phase. Even smaller halos never reach the threshold 
mass for quasar triggering; these still contain star-forming disk galaxies at z < 0.8, and we observe some of them as optical or IR-selected Seyfert galaxies. The 
dashed box indicates the AGN types (in their characteristic dark matter halos) that would be observable in the redshift range 0.25 < z < 0.8. 



mergers (with m ass ratio s >3:1) wi th other gas rich systems 
dRobertson et al . 2006; Ho pkins et a l. 2009a). As the galaxies 
grow, their dark matter halos also increase in mass through 
the hierarchical growth of structure. The largest halos grow 
fastest, so that the most massive overdensities at high red- 
shift remain the most massive at lower redshift, a nd so on for 
smaller masses (e.g.. Press & Schechter 1974; La cey & Cold 
[1993). 

When a galaxy's parent dark matter halo reaches a criti- 
cal Mhaio between 10'^ and lO'"^ M©, a dramatic event occurs 
that triggers luminous quasar activity and thus rapid growth 
of the central black hole, as well as producing a dynami- 
cally hot stellar system (the classical bulge). One candidate 
for such an event is the major merger of two gas rich galax- 
ies; the characteristic quasar halo mass of Mhaio 3 x 10'^ 
Mq corresponds to the small group environment where ma- 
jor m ergers of galaxies with masses ~ are most com- 
mon (' Hopkins et aT] l2008b). Models of mergers appear to 
naturally produce many observed properties of quasars and 
the resulting spheroidal remnants, as well as their abun- 



dance and r edshift evolution ("Kauffmann & Haehneltl 120001: 
ISpringel eta l. 2005; Hopkins et al. 2008b), and there is evi- 
dence for merger remnants hos ting individual quasars (e.g., 
iRiechers et allboOBl: lAravena et aL.2008;) . 

However, in addition to mergers, other typ es of events such 
as disk instabilities (e.g., i Bower et al.ll2006l) may also possi- 
bly trigger efficient accretion onto the black hole and grow the 
stellar bulge. We also note that other processe s, such as the ac- 
cretio n of recycled gas from evolved stars ( Ciotti & Ostrikerl 
l2007h may trigger a quasar after the formation of the bulge. 
In the context of our observations at z ~ 0.5, the details of 
the triggering mechanism are relatively unimportant. It is suf- 
ficient only that the primary mechanism for quasar accretion 
and the growth of the stellar bulge occur together at a charac- 
teristic dark matter halo mass, at high redshift for the largest 
overdensities, and at progressively later times for smaller and 
smaller overdensities. 

After the event that builds the bulge and grows the central 
black hole, star formation in the galaxy must be quenched 
in order to produce the observed population of passively 
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evolving, bulge-dominated galaxies with old ste llar popula- 
tions that are found on th e red sequence (e.g., iNelan et aP 

f\03; 'Thomas et all l2005h. \yhich i s detected out to z ~ 2 
;irmetal. 20081: iKriek et al.ll2008l) . The cessation of star 
formation must necessarily occur after or simultaneous to the 
formation of the bulge, since essent ially all pa ssively evolv- 
ing galaxies contain stellar bulges ( lBellll2008l) ^^. It is pos- 
sible that this quenching is a direct result of the event that 
fuels the quasar. Processes that drive gas onto the black 
hole can also power intense star bur sts (as are often as- 
sociated with lumi nous AGNs, e.g. Alexander et al.l 120051; 
iBrand et al.1 r2006bl) that can use up the galaxy's supply of 
cold gas and thus halt star formation. In addition, high 
virial temperature of the (growing) dark matter halo can 
limit further accretion of cold gas and thus star formation 
(e.g. Rees & Ostriker 1977; Cattaneo et al. 2006; Fa ber et all 
l2007h . Simulations suggest that the tra nsition to this "hot 
halo" regime occurs at Mhaio ~ 10'^ (e.g.. lBirnboim & Dekell 
l2003HDekel & Birnb oim 2006), close to the typical Mhaio for 
quasar activity. For our simple picture, the details of this 
quenching are not important; it is sufficient that soon after 
the formation of a massive stellar bulge, star formation in the 
galaxy rapidly ceases. 

If the timescale for the quenching of star formation is short, 
the galaxy's optical spectrum will be a composite of emis- 
sion from the old stellar population plus a fading contribu- 
tion from the younger stars, producing "green" colors that 
become redder as the young stellar population evolves. The 
timescale for the galaxy to evolve red colors typical of early- 
type galaxies is ~l -2Gy r (iNewberry et al.ll990tlBarger et al.l 
[1996; Bow er etal.H l998|). This phase may be accompanied 
by further accretion onto the central black hole, which can 
evolve on a similar timescale (iHopkins et al.ll2005l) . and may 
contribute to the galaxy's transition to the red sequence. Ob- 
servations of local (z < 0.1) early-type galaxies suggest that 
optically-detected AGN activity coincides with the decline 
of star formation and th e transition from blue to red colors 
dSchawinski et al.'l2 007[). as well as t he destruction of molec- 
ular clouds (Schawii iski et al.ll2009al) . The decline of the av- 
erage black hole accretion rate from its peak in the quasar 
phase (when M ~ MEdd) can occu r (with significant y ari- 
ability) over timescales of ^^1 Gyr ( iHopkins et al.ll2005h . fi- 
nally reaching accretion rates of M < lO'^Mndd or lower At 
these low accretion rates (and corresponding low luminosi- 
ties), the S EP of an A GN is increasingly dominated by X-ray 
emission (Steffen et al. 2006 ; Vasudeva n & Fabian 2007), and 
may even enter a new, radi atively inefficient mode of accre- 
tion tehurazov et al.ll2005h . Therefore, these fading AGNs 
in "green" galaxies may be most efficiently found as X-ray 
AGNs. 

AGN activity may serve to prevent further bursts of star for- 
mation. As mentioned above, at this stage the massive stellar 
bulge of the galaxy will be surrounded by a halo of hot, viri- 
alized gas. While most of this gas will have very long cool- 
ing times, near the center (where the gas density is highest 
and cooling most efficient), gas cooling could fuel further star 
formation. Such cooling must be suppressed to explain the 
old stellar populations observed in red-sequence galaxies and 
the obs erved dearth of m assive, blue star-forming galaxies at 
z < 1 (Fa ber et aP l2007h . Once cooling is suppressed, the 
red-sequence galaxies can further grow through "dry" merg- 

We note that the lBelll j20O8h result is only for central galaxies in groups 
or clusters; see Section l9i4l f or a discussion. 



ers that would not induce further bursts of star formation. 

One promising mechanism for suppressing cooling in 
red sequence galaxies is feedback from AGNs, which 
can inject mechanical e n ergy t hat reheats the cool- 
ing gas (Tabor & Bhine^ 119931: iBinney & Tabo3 119951: 
iChurazov et al. 2002). While there is indirect evidence for 
such feedback from op t ically or X-ray-sele cted AGNs (e.g., 
ISchawinski et al.1 120071: iBundv et al.1 l2008l) . the process of 
AGN feedback is observed most directly in radio AGNs. 
In many elliptical galaxies and galaxy clusters, relativis- 
tic jets from AGNs (detected in the radio) inflate buoyant 
bubble s and drive shocks in the hot surrounding medium 
(e.g i Fabian et al.'2003VMcN ainara et al.l2005HForman et all 
l2007M McNamara & Nulsen 20Q7|)- The heating rate from this 
process can be sufficient to balance the cooling from t he hot 
halo, and might t herefore q uench star formatio n ( BirzaiT et al.l 
l2004HRaffertvet al. 2006; ICroston et al.ll2008l) . The Edding- 
ton ratios for accretion in these "radio-mode" outbursts are 
generally small (A < 10"^), and the observed radiative power 
is far smaller than the inferred mechanical power of the jets 
(Churazov et al. 2002; Birzan et al. 2004). It has been pro- 
posed that these radiatively inefficient, jet-dominated out- 
bursts may be fueled by accretion directly from the hot gas 
halo (in contrast to the accretion of cold gas that likely pow- 
ers high-Eddington, optically luminous AGNs), and so is only 
possible in massive galaxies with large dark matter halos (e.g., 
Tasse 2008). 

In summary, the simple picture described above suggests 
that massive galaxies undergo the following evolutionary se- 
quence: gas-rich, star-forming, rotation-dominated system 
SMG/ULIRG/quasar — + "green" spheroid with declining 
AGN activity — > re d spheroid with intermittent radio AGN ac- 
tivity. Figure 1 of IHopkins et al.l (l2008b ) shows a schematic 
of these various phases, for the case in which the quasar is 
triggered by major mergers. The SMG/ULIRG/quasar phase 
evolves on a short timescale (<10^ yr) and occurs in dark mat- 
ter halos with a characteristic Mhaio between 10'^ and lO'-' 
Mq. The subsequent "green" spheroid phase evolves over 
somewhat longer times (~1 Gyr), while the phase of inter- 
mittent radio AGN activity can last for the subsequent life- 
time of the system. A schematic of this sequence is given 
in Figure [16] The columns represent systems for different 
initial halo masses, for which the evolutionary sequence oc- 
curs at different redshifts: large dark matter halos reach the 
critical mass lO'^-lO'-' Mq (and thus trigger quasars) at rel- 
atively high redshifts, while for smaller systems the sequence 
occurs later Therefore, at z ^ 0.5, objects in different stages 
of this sequence should be in dark matter halos of different 
sizes, which can be probed through their clustering. 

A key element of this simple model is that quasar activity, 
bulge formation, and the quenching of star formation occur 
when dark matter halos reach the critical Mhaio- Therefore, 
galaxies residing in halos that have not reached this mass will 
remain as gas-rich disks, and will not experience a luminous 
quasar phase and the associated rapid black hole growth. A 
schematic of the evolution of low halo mass systems is given 
in the rightmost column of Figure [16] Although these objects 
have never hosted luminous quasars, they can be seen as low- 
luminosity AGNs; various processes such as the stoch astic 
accretion of cold gas clouds (IHopk ins & Hernquistll2006 ) can 
fuel low-level (and likely intermittent) AGN activity, while 
galaxy-scale processes, such as the formation of galactic bars, 
can enhance the gas density in the nuclear regions and prefer- 
entially help drive accretion (e.g.. lOhta et al.ll2007h . 
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Accretion in low halo mass systems can proceed in a wide 
range of accretion rates; however, owing to the small black 
hole masses (Mbh S ^0) only the sources with rela- 
tively high accretion rates will be detected at z 0.5 in a 
wide-field survey such as Bootes. The observational evidence 
for AGN "downsizing", in which small black holes domi- 
nate the accretion density of the univers e at low redshift (e.g., 
iHasinger et"ani2005l; iBarger et al.ll2005h . suggests that AGNs 
in less-luminous, disk-dominated galaxies should dominate 
the population of high-Eddington-rate AGNs at z < 0.8. These 
high-Eddington AGNs would likely be identifiable by their 
mid-IR colors. 

9.3. Comparison to observations 

We next compare the evolutionary scenario described above 
to our observational results for the host galaxies, clustering, 
and Eddington ratios of the different classes of AGNs. 

1 . Radio AGNs. Most radio AGNs in our survey at 0.25 < 
z < 0.8 are found in luminous red sequence galaxies 
and are strongly clustered. In the evolutionary picture 
described above, radio AGNs are a key part of the late 
stages of massive galaxy evolution, after the bulge has 
formed and star formation has ceased. Therefore at 
z ^ 0.5, radio AGNs would be expected to be found pri- 
marily in massive red sequence galaxies in large dark 
matter halos. Accordingly, we find that radio AGNs 
are located in relatively large dark matter halos with 
Mhnio ~ 3 X lO'"* Mq, corresponding to the large 
galaxy groups or small clusters. The radio AGNs are 
clustered similarly to a control sample of galaxies with 
matched luminosities and colors, suggesting that inter- 
mittent radio activity may be common, at least for the 
host galaxies and radio luminosities probed by our sam- 
ple. 

We note that a small fraction 20%) of the radio 
AGNs with extended optical counterparts are found in 
blue galaxies, which is contrary to this simple pic- 
ture. However, some of these are are likely radio- 
loud Seyfert galaxies and quasars, which make up 
~10% of those populations, and which may represent a 
mode of accretion more similar to radiatively-efficient 
Seyf erts than mechanically-dominated radio ga laxies 
(e.g.. [Merloni & Hein3l2008l: [Ponoso et al.ll2009l) . 

2. X-ray AGNs. X-ray AGNs in our survey at 0.25 < z < 
0.8 are found in host galaxies throughout the color- 
magnitude diagram, but the distribution peaks in the 
"green valley" along with a tail of bluer galaxies. The 
large-scale (1-10 Mpc) clustering indicates that the 
X-ray AGNs inhabit dark matter halos of mass ^ lO'-' 

Mq, similar to the clustering for normal AGES 
galaxies. The X-ray AGNs are unbiased relative to 
a control sample with matched galaxy properties, al- 
though there is a (marginally-significant) antibias rel- 
ative to AGES galaxies on scales < 1 /i"' Mpc, sug- 
gesting that the X-ray AGNs may preferentially reside 
in central galaxies. In our simple evolutionary picture, 
central galaxies in halos with Mhaio ~ lO'-' /i"' Mq are 
likely to have recently experienced the buildup of the 
stellar bulge and quenching of star formation, followed 
by a decline in the accretion rate onto the black hole. 
The properties of X-ray AGNs, with hosts typically 



showing "green" colors and with a wide range of Ed- 
dington ratios (10~^ < A < 1) are generally consistent 
with this picture. 

3. IR AGNs. IR-selected AGNs in our observations are 
found in relatively low-luminosity galaxies and are 
weakly clustered, with characteristic Mhaio ^ 10'^ 
Mq. There is a significant antibias relative to a matched 
control sample of normal galaxies. The Eddington ra- 
tios of these sources are relatively high (A > 10"^), 
and the clustering and Eddington ratios do not depend 
strongly on whether the AGNs are also detected in 
X-rays. In our simple picture of AGN evolution, at 
z < 0.8 most high-Eddington AGNs will be found in ob- 
jects with small black holes, in environments that have 
not yet reached the critical halo mass for the quench- 
ing of star formation. In our observations we identify 
IR-selected AGNs with weakly clustered star-forming 
galaxies, in general agreement with this picture. 

9.4. Caveats 

The simple evolutionary scenario described above thus cor- 
responds generally to our observational results of radio, X- 
ray, and IR-selected AGNs. However, by its nature our simple 
model of AGN and galaxy evolution ignores many details and 
is subject to several important caveats. In this section we dis- 
cuss several of these issues and how they may be addressed 
with future observations. 

1 . Central versus satellite galaxies. The model described 
above may be only strictly valid for central galaxies, 
defined as the dominant galaxy in the parent dark mat- 
ter halo. It is possible that the quenching of star for- 
mation can happen in different ways for satellite galax- 
ies that fall into larger dark matter halos; ram pressure 
stripping, strangulation, and starvation can cause a disk 
galaxy to lose its cold gas and cease star formation (see 
Ivan den Bosch et alJ2008 . and references therein). Fur- 
ther, there is evidence that for satellite galaxies, star 
formation and AGN activity depend more strongly on 
the stellar mass of th e galaxy than the mass of the halo 
jPasquali et al.ll2009l) . Successful models of AGN and 
galaxy evolution will likely need to explain the evo- 
lution of central and satellite galaxies separately. We 
note, however, that at least for X-ray and IR-selected 
AGN in our sample, the observed clustering suggests 
that they may preferentially reside in central galaxies 
(SectionlT2]i. 

2. Critical halo mass. Our simple evolutionary model in- 
cludes "critical" dark matter halo mass Merit between 
10'^ and 10'^ Mq at which a galaxy undergoes lu- 
minous quasar activity and the growth of the stellar 
bulge. This idea is motivated by observations of quasar 
clustering and models of hot gaseous halos, but the 
details of this process are certainly more complicated 
than can be described by a single critical halo mass. 
Indeed, galaxy evolution models in which quench- 
ing occurs at a halo mass threshold may fail to re- 
produce the density of passively-evolving galaxies at 
high redshift, as well as relative numbers of red and 
blue galaxies ( Hopkins e t al. 2008a). Further, evidence 
from the evolution of luminosity function on the red 
and blue sequence suggests that quenching may occur 
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with different mec hanisms, and in different mass halos 
(iFaber et al.ll2007h . Quasar triggering and bulge forma- 
tion probably occur by some more complex process that 
preferentially occurs at a given halo mass. The details 
of this process cannot be constrained with the obser- 
vations presented here, but they could be tested in the 
future by more detailed work. 

3. Extinction and galaxy colors. The above discussion 
makes note of the "green" galaxy host colors of the 
X-ray AGNs as evidence that these sources have un- 
dergone recent star formation. However, in calcu- 
lation of galaxy colors, we have not taken into ac- 
count the effects of dust extinction intrinsic to the 
galaxy. It is possible that some galaxies in the "green 
valley" are in fact actively star-forming galaxies that 
are reddened by dust. In a multiwavelength study of 
galaxies at 0.05 < z < 1.5 in the GOODS-North field, 
ICowie & Barged (l2008h found a significant population 
of galaxies (mainly Spitzer 24 /im sources) that have 
observed colors in the "green valley", but lie in the blue 
cloud after correcting for extinction. However, in our 
sample the X-ray AGNs in "green" galaxies show clus- 
tering consistent with that of typical quiescent "green" 
galaxies with similar colors (Section l7T2] ). and stronger 
clustering than typical blue cloud galaxies. Thus it is 
unlikely that the majority of the "green" hosts of X-ray 
AGNs are dust-reddened blue cloud galaxies. In the 
future, a more detailed study of the optical-IR SEDs 
and average optical spectra of these objects would al- 
low us to put more strict limits on the population of 
dust-reddened star-forming galaxies. 

10. SUMMARY 

In this paper we explored the links between SMBH ac- 
cretion, host galaxy evolution, and large-scale environment, 
by studying the host galaxies of AGNs detected at redshifts 
0.25 < z < 0.8 in the AGES survey. The AGNs were selected 
using observations in the radio. X-ray, and IR wavebands. 
Key results from this paper are as follows: 

1. AGNs in the AGES sample at 0.25 < z < 0.8 that are 
selected in different wavebands comprise generally dis- 
tinct populations of sources. Radio AGNs are generally 
not selected in the other bands, while there is 30%-50% 
overlap between the X-ray and IR-selected AGNs. 

2. The host galaxies of optically faint AGNs are signifi- 
cantly different between the three AGN samples. Radio 
AGNs are found in luminous red-sequence galaxies. X- 
ray AGNs are found in galaxies of all colors, with a 
peak in the "green valley". IR AGN hosts are relatively 
bluer and less luminous than those of the X-ray or radio 
AGNs. 

3. The two-point cross-correlation between AGNs and 
galaxies is significantly different for the three classes of 
AGNs. We estimate that radio. X-ray, and IR AGNs are 
found in dark matter halos with characteristic masses 
-3 X 10'^ /i"' Mq, ~10'^ /i-i Mq, and <10'2 Mq, 
respectively. 

4. X-ray and radio AGNs with extended galaxy coun- 
terparts are clustered similarly to samples matched in 



color, absolute magnitude, and redshift, indicating that 
they inhabit environments typical of their parent host 
galaxy populations. Conversely, IR AGNs are weakly 
clustered relative to a matched galaxy sample, suggest- 
ing that the process that triggers IR-bright AGN ac- 
cretion depends on local environment as well as host 
galaxy type, and is more hkely to occur in regions of 
lower galaxy density. 

5. The average X-ray spectra of X-ray and IR AGNs are 
consistent with the properties of their optical counter- 
parts. Optically bright AGNs have soft X-ray spec- 
tra with photon index F ~ 1.6, typical of unabsorbed 
AGNs. Optically faint AGNs have harder average spec- 
tra, which likely indicates either absorption by interven- 
ing gas, or a radiatively inefficient mode of accretion. 

6. Most radio AGNs in our sample have massive black 
holes (Mbh > 10^ Mq) and very small Eddington ra- 
tios (A < 10"-'). X-ray AGNs have smaller typical Mbh 
values that extend down to ^^10^ Mq, and A between 
10"^^ and 1. IR AGNs have relatively small black holes 
(3 X 10^ < Mbh < 3 x 10*^ Mq) and high Eddington ra- 
tios (A > 10-2). 

7. Sources that are selected as both X-ray and IR AGNs 
have host galaxy properties, clustering, and Eddington 
ratios that are similar to those for the full IR AGN pop- 
ulation. This suggests that the IR-bright AGN may rep- 
resent a single population of AGNs, distinct from typ- 
ical X-ray AGNs (most of which are not selected with 
IRAC). 

8. The host galaxies, clustering, and Eddington ratios of 
the three classes of AGNs are generally consistent with 
a simple picture of AGN and galaxy evolution in which 
SMBH growth and the quenching of star formation oc- 
cur preferentially at a characteristic dark matter halo 
mass Mhaio ^ lO'^-lO'-' Mq, and are followed by a de- 
cline in the AGN accretion rate and a change to a ra- 
diatively inefficient, mechanically dominated accretion 
mode. 
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Fig. 17. — Average flux profiles for subsets of galaxies and AGNs from the AGES survey, in the Bw (top) and / (bottom) bands. Fluxes are calculated in 
annular apertures of width 1" (the A-axis gives the outer diameter of the annulus). Shown are the profiles for normal AGES galaxies (black solid line), AGNs 
with extended optical counterparts (radio, orange dotted fine; X-ray, green solid line; IR, red dashed line). X-ray AGNs with unresolved optical counterparts 
(blue dot-dashed line), and a 1"35 (FWHM) Gaussian PSF (black dotted line). The AGNs in extended sources have relatively flat flux profiles simUar to normal 
galaxies, suggesting that they have little nuclear contamination. 



APPENDIX 

A. ESTIMATING NUCLEAR CONTAMINATION FOR AGN HOST GALAXIES 

Here we describe in detail our method to estimate the contribution of active nuclei to the optical colors of X-ray and radio- 
selected AGNs. For unobscured, optically bright Seyfert galaxies and quasars, the optical emission from the nucleus is generally 
significantly bluer than the host galaxy, so that the combined AGN plus galaxy is bluer than the galaxy alone. For our study 
of the colors and luminosities of AGN host galaxies (Figure |9]l we have chosen sources that are optically extended (and so not 
dominated by a central point source), but the nucleus may make some contribution to their integrated optical flux. 

To test for nuclear contamination, we use the optical imaging to calculate the optical flux and color profiles of the different 
types of sources. Here we use aperture photometry derived using the technique of Brown et al. (2007, 2008), with the Bw, R, and 
/-band photometry smoothed to a common PSF with FWHM 1 ."35 (to avoid PSF differences between bands), and uncertainties 
estimated using a Monte Carlo technique. 

Evidence for some nuclear contamination is given in Figure [T71 which shows the average Bw and /-band flux profiles for 
normal galaxies, and radio. X-ray, and IR AGNs with extended optical counterparts. For comparison, we also show profiles for 
X-ray sources with unresolved optical morphologies (see Section |4!4] |. as well as the expected profile for a Gaussian PSF with 
FWHM 1"35, typical of the NDWFS observations. The profiles are calculated in annuli of 1" width. Normal galaxies have 
broad, extended profiles, while optically unresolved X-ray AGNs have more strongly peaked profiles indicating that they have a 
significant contribution from a nuclear point source. The radio AGNs have a broad average profile similar to that of quiescent 
galaxies, indicating that on average, they have little or no nuclear contamination (this may be expected, given that most radio 
AGNs also show little nuclear activity in the X-rays or IR). In contrast. X-ray and IR-selected AGNs show a slight excess in the 
Bw band over the quiescent galaxies inside a 3" diameter, suggesting a moderate contribution from a nuclear point source. There 
is only a weak excess for the X-ray and IR AGNs in the / band, however, indicating that the nuclear contribution is bluer than the 
host galaxy. Therefore nuclear emission may affect the integrated colors of the source. 

For the radio. X-ray, and IR AGNs, we wish to estimate the extent of this optical color contamination by the active nucleus, 
and thus produce corrected rest-frame '(w-r) colors for the host galaxies alone. To this end, we compare the optical colors 
observed in the center of AGN-hosting galaxies to that at the outskirts, and compare these to quiescent galaxies. First, for all 
AGES galaxies we derive the quantity C3-7, which is the Bw-I color in an annular aperture of inner and outer diameters 3" and 
7", respectively (Figure [Tsll. Since the observations have a PSF of width L"35 (FWHM), excluding the central 3" removes most 
of the nuclear contribution while including enough flux from the outskirts of the galaxy to derive sufficiently accurate colors. We 
then calculate AC, the difference between Ca-v and C4 (the Bw-I 4" aperture colors). We would expect nuclear contamination 
to make C4 bluer, but to have little effect on Cj,--i, so that contamination should increase AC. 

We note that the color profiles of galaxies vary with galaxy type; early-type galaxies are redder in the centers (where they 
are dominated by old, higher metallicity stars) compared to late-type galaxies. Figure [19] shows AC versus the rest-frame 
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Fig. 18. — NDWFS /-band image of a randomly selected X-ray AGN with an extended optical counterpart at z = 0.55. The circles show the apertures used for 
determining the color contamination from the nucleus; the solid blue line shows the inner 4" aperture, while the red dashed lines show the annulus with inner 
and outer radii of 3" and 7", respectively. 




1 2 3 1 2 3 

° '(M - r) - r) 

Fig. 19. — Color profiles of galaxies and AGNs in the AGES sample. Shown is AC (the difference in Bw —I color between the 3"-7" annulus and the central 
4" aperture) vs. rest-frame color " '(w— r), as derived in Section [33l The contours and points represent normal galaxies, while the open circles, stars, and squares 
represent radio, X-ray, and IR AGNs with extended optical counterparts. The bottom right panel shows the same quantities for AGNs with broad lines in their 
optical spectra. Optically unresolved sources are shown in dark blue, and optically extended sources are shown in fight blue. In all panels, the thick lines sfiow 
linear fits to AC vs. " '(« — /■) for the three sets of sources (see the text for details). 

" '(M-r); contours and points show this distribution for normal galaxies, while stars and circles represent X-ray and radio AGNs, 
respectively. For normal galaxies, AC increases slightly with bluer " '(M-r); the reddest galaxies have AC ~ -0.1 while the 
bluest have AC ~ 0. In contrast, the radio. X-ray, and IRAC AGNs show stronger increases in AC with bluer colors, indicating 
that some X-ray AGNs have significant nuclear emission. The solid black and dashed green, orange, and red lines show linear 
average fits to the trend of AC with " '(m-t). The trends are roughly linear; there is no significant difference if we compute the 
average AC in bins of " '(M-r), so for simplicity we use the linear fit. 

Using these fits, we can derive the average nuclear contamination of the observed color as a function of the integrated " '(M-r) 
obtained for X-ray and radio AGN. For a given " '(M-r), the Byi/-I A" color of the host galaxy (Cgai) is related to the observed 
4" color (C4) by 
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Fig. 20. — A'-corrected rest-frame color (w — r) v.s. observed Bw -I 4" aperture color, for AGES main sample galaxies (points), and AGNs with extended 
optical counterparts: radio (circles). X-ray (stars), and IRAC (squares). Sources and broken power-law fits are shown for two redshift bins. 



Cgal=C4 + (ACW, (Al) 

where (AC)corr is a function of " '(m- r), and is simply the difference between the dashed green line and the solid black line in 
Figure[T9] Including all radio. X-ray, and IR AGNs, the best linear fit gives: 

( AQcon- = 0.58 + 0.20['''(m - r)] . (A2) 

We note that, as is clear in Figure[T9] (AC)coit versus " '(M-r) is slightly different for the three classes of AGNs, but for simplicity 
we use the fit that includes for all AGNs. 

For normal galaxies and AGNs, the average AC increases slightly with redshift; if we divide the sample into sources with 
0.25 < z < 0.5 and 0.5 < z < 0.8, the average AC differs by up to ~0.2 (depending on " '(w-r)). However, the change with 
redshift is the same for both normal galaxies and AGNs, so (AC)(;oit [as a function of ° '(m - r)] changes by very little (less than 
0.04). For simplicity, we therefore use the single parameterization for (AC)coit versus " '(M-r) (Eqn. lA2l i for all 0.25 < z < 0.8. 

For comparison to the radio. X-ray, and IR AGNs with extended counterparts, we also derive AC for X-ray sources at 0.25 < 
z < 0.8 with broad emission lines in their AGES optical spectra, indicating that they are optically unobscured AGNs'^. These 
objects are shown in the lower right panel of Figure [19] Sources with unresolved morphologies (broad-line, unresolved AGNs) 
are shown in dark blue, and those with extended morphologies are shown in light blue. The broad-line, optically unresolved 
AGNs typically have AC ~ 0.5 and have " '(m- r) < 1, and they do not show a trend in AC with " '(m- r), while the optically 
extended sources show a similar trend to the radio. X-ray, and IR AGNs in the other four panels, although the broad-line AGNs 
have slightly stronger average contamination. Since most AGNs with extended galaxy counterparts have smaller AC and redder 

'(M-r) than is typical of optically unresolved AGNs with broad lines, we conclude that these sources are not as dominated by a 
nucleus (as also indicated in Figure[T7l and that their colors may be robustly corrected for nuclear contamination. 

To correct the colors for nuclear contamination, we use Eqn. lA2l to derive (AC)coit for each source, and so determine a "cor- 
rected" C4, then convert this to a corrected rest-frame color " '(m-t). To determine how " '(M-r) varies with C4, we divide AGES 
main sample galaxies into seven bins of redshift from 0.25 to 0.8, and fit the dependence of °'(M-r) on C4 in each bin. A broken 
power-law model with a break at C4 = 3.3 gives a good fit for all redshifts; examples of two fits are given in Figure |20] The 
relation between " '(«- r) and C4 becomes shallower with increasing redshift, with slopes given by: 

A" '(M-r) _ [0.73+ 1.27(z- 0.70)^ (C4 < 3.3) 

AC4 \0.33 + 3. 17(z- 0.72)2 (C4 > 3.3) ^ ' 

We thus calculate the correction to '(M-r) by multiplying (AC)con (from Eqn. IA2l i by A"^(M-r)/AC4 (from Eqn. lA3l l. The 
typical corrections are small, varying from ~0.3 for the bluest AGNs, to less than 0.05 for the reddest AGNs. The effects of these 
nuclear contamination corrections on the host-galaxy color distribution of the AGN sample are shown in Figure|9l 

As mentioned in Section [T] in this paper we generally do not make use of optical spectroscopic data for selecting AGNs, because of significant selection 
effects. In estimating nuclear color contamination, however, these sources are useful because the spectroscopy provides evidence for some optical contribution 
from the AGN. 
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log(Source counts [0.5-7 keV]) Source counts (0.5-7 keV) 

Fig. 21 . — (a) Histogram of source counts for the 50 IR AGNs with extended optical counteipails that ai'e detected in X-rays. We extrapolate the roughly 
log-uniform disMbution to lower fluxes, to model the observed counts for the X-ray undetected sources, (b) Histogram of total observed counts (source plus 
background) within rjo for the 82 IR AGNs with extended optical counterparts that are inside the XBootes coverage area and are not detected in X-rays. The red 
line shows the predicted counts distribution from a simple extrapolation of a log-uniform flux distribution observed for the detected sources. 



B. LIMITS ON STAR BURST CONTAMINATION OF IR-SELECTED AGN SAMPLE 

As discussed in Section l431 it is possible that the IRAC-selected AGN sample may contain some sources whose IR emission is 
powered by star formation instead of nuclear accretion. For the IRAC AGNs that are also detected in X-rays (roughly half of the 
total sample, and^l /3 of those with extended optical counterparts), the X-ray emission is an unambiguous indicator of nuclear 
accretion. In addition, ^40% of IRAC-selected AGNs show clear evidence of nuclear emission on the basis of their unresolved 
optical morphologies. However, for 83 (or 30%) of the IR-selected AGNs that are optically extended and not detected in X-rays, 
it is useful to estimate the possible fraction of star burst contaminants. 

For this we also u se X -ray emission, but are limited to X-ray stacking because the sources are not individually detected. As 
discussed in Section [STl the X-ray undetected IRAC AGNs have a high average luminosity i{Lx) ~ lO'*^ erg s"') and hard X- 
ray spectrum with HR = 0.21 ±0.21 (or F ~ 0.2), characteristic of absorbed AGNs. The stacked signal is significantly harder 
and more lu minous than the typical X-ray spectra of star burst galaxies, which have Lx ^3 x 10*^ erg s"' and F > 1.4 (e.g., 
iRanalli eral.ll2003l) . This indicates that most of the stacked X-ray flux from these sources is from AGNs, but there remains a 
possibility than an X-ray-faint population of star bursts contributes significantly to the stacked flux. 

To test this possibility, we extrapolate the X-ray flux distribution for the IRAC AGNs with extended counterparts that are X-ray 
detected, and check that the stacked signal is consistent with this extrapolation. As a constraint, we use the distribution of source 
counts at each IRAC AGN position. This analysis constrains not only the ave rage X-ray flux of a sample, but also the shape of 
its underlying flux distribution (for a detailed discussion of this technique, see Hickox & Markevitchll200'7l) . 

The X-ray flux distribution for the 50 IRAC-selected AGNs with extended counterparts that are detected in X-rays is roughly 
log-uniform (Figurel2Tl;). The IR luminosity and redshift distributions of the X-ray- detected sources are similar to those of their 
X-ray undetected counterparts, so it is reasonable to expect that the shape of the X-ray flux distribution of the IRAC AGNs may 
extend to fainter fluxes. We therefore create a log-uniform model flux distribution for the 82 X-ray undetected IRAC AGNs that 
have extended optical counterparts and lie within the X-ray coverage area. The model sources have fluxes from 0.02 to 4 count 
source"' (29 sources dex"' in flux). We randomly assign the model fluxes to each source in the sample. 

We next calculate an estimated background X-ray flux from each undetected IRAC AGN, given the average background surface 
brightness in the 0.5-7 keV band (Section lOT l. and the area of the extraction aperture (given by the 90% energy encircled radius). 
The background varies from 0.01 to 0.1 count source"'. Adding the background to the flux from the model flux distribution, we 
calculate a predicted distribution of source counts at the IRAC AGN positions. We repeat this calculation 10 times, with different 
permutations of the source and background fluxes, and average the results to produce a model counts distribution. 

This model distribution shows remarkably good agreement with the total number of observed source counts and their distribu- 
tion, as shown in Figurel^Tb (note that the model is a direct extrapolation of the observed flux distribution, and has not been scaled 
fit the observed counts). The X-ray emission from the undetected IRAC AGNs is therefore fully consistent with an extrapolation 
of the observed flux distribution for the X-ray-detected sources. 

Of course, there remains the possibility of a population of star-forming galaxies at the faint end of the sample. To put limits 
on such a population, we repeat the above calculation but replace the faintest model AGNs with star burst galaxies, with X-ray 
fluxes of 0.05 c ount source"' in 5 ks (corres ponding to Lx ^ 10"" erg s"'). We calculate the variation in the C-statistic llCashI 
(Il979h . see also lHickox & Markevitclil (I2007M as a function of the fraction of sources replaced by star bursts. The best match to 
the data is achieved for the case of no star burst contamination; to 95% confidence, the non-AGN fraction must be less than 60%. 
This implies that for the full sample of IRAC AGNs with extended optical counterparts, the contamination is less than 40%. 

In addition to the X-ray constraints on contamination, we also note that (1) these sources would have to be highly extinguished 
star bursts in order to have the AGN-like IRAC colors, and (2) the optically extended IR AGNs show small but significant 
nuclear contamination in the optical. On the whole, these facts suggest that the star burst contamination of the IRAC AGNs 
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with extended counterparts is most likely quite small (<20%). This implies that among the full IRAC AGN sample (including 
optically unresolved sources), star burst contamination is at most 20% and more likely <10%. 
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